A Class of P, T-Invariant Topological Phases of Interacting Electrons 
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We describe a class of parity- and time-reversal-invariant topological states of matter which can 
arise in correlated electron systems in 2+1-dimensions. These states are characterized by particle-like 
excitations exhibiting exotic braiding statistics. P and T invariance are maintained by a 'doubling' 
of the low-energy degrees of freedom which occurs naturally without doubling the underlying mi- 
croscopic degrees of freedom. The simplest examples have been the subject of considerable interest 
as proposed mechanisms for high-T c superconductivity. One is the 'doubled' version of the chiral 
spin liquid. The chiral spin liquid gives rise to anyon superconductivity at finite doping and the 
corresponding field theory is (7(1) Chern-Simons theory at coupling constant m = 2. The 'doubled' 
theory is two copies of this theory, one with m = 2 the other with m = —2. The second example 
corresponds to Z2 gauge theory, which describes a scenario for spin-charge separation. Our main 
concern, with an eye towards applications to quantum computation, are richer models which support 
non-Abelian statistics. All of these models, richer or poorer, lie in a tightly-organized discrete family 
indexed by the Baraha numbers, 2 cos j^pg, for positive integer k. The physical inference is that a 
material manifesting the Z2 gauge theory or a doubled chiral spin liquid might be easily altered 
to one capable of universal quantum computation. These phases of matter have a field-theoretic 
description in terms of gauge theories which, in their infrared limits, are topological field theories. 
We motivate these gauge theories using a parton model or slave-fermion construction and show 
how they can be solved exactly. The structure of the resulting Hilbert spaces can be understood 
in purely combinatorial terms. The highly-constrained nature of this combinatorial construction, 
phrased in the language of the topology of curves on surfaces, lays the groundwork for a strategy 
for constructing microscopic lattice models which give rise to these phases. 



I. INTRODUCTION 

In quantum mechanics, trajectories in which identical 
particles are permuted must be considered on an equal 
footing with those in which they are not. Consequently, 
particles can be classified according to the irreducible 
representations of the permutation group. This leaves 
bosons and fermions as the only allowed possibilities 1 . 
However, in two spatial dimensions, trajectories can be 
further classified according to the braid group, essentially 
because the combination of two counter-clockwise ex- 
changes in succession can not be adiabatically deformed 
to no exchange at all. As a result, non-trivial braiding 
statistics is possible in two dimensions^. 

The simplest examples of particles with exotic braid- 
ing statistics are called anyons. They realize one- 
dimensional representations of the braid group: the 
quantum-mechanical wavefunction acquires a non-trivial 
phase as a result of a counter-clockwise exchange of one 
anyon with another (the complex conjugate phase is as- 
sociated with a clockwise rotation) or the 2tt rotation of 
an anyon. It is actually not necessary for the anyons 
to be identical, unlike with the permutation group. For 
instance, a phase e % ®" could result when a particle of 
type i is exchanged with another of type i while a phase 
e y results when a particle of type i winds around a 
particle of type j and returns to its original position 
(when the particles are distinguishable, a trajectory in 
which they are exchanged leads to a different configu- 
ration). Thus, the term 'statistics' is somewhat mislead- 



ing in the two-dimensional case because the classification 
is not according to the permutation group. It is more 
correct to say, instead, that there is a topological inter- 
action at work: an interaction between particles which 
only depends on how they are braided. This is reflected 
in the field-theoretic implementation of anyonic statis- 
tics, in which a (7(1) gauge field mediates the interac- 
tion between the particles. If the gauge field is governed 
by the U(l) Chern-Simons action, then the interaction 
is purely topological in nature. The coefficient of the 
Chern-Simons term determines the phase assigned to a 
particle exchange. The quasiparticlc excitations of most 
observed fractional quantum Hall states are believed to 
be anyons^i^. 

An even more exotic variety of braiding statistics is 
associated with multi-dimensional representations of the 
braid group. If there are p states of the system when 
there are particles at Xi, X2, . . . , x„, then the effect of an 
exchange of identical particles may be represented by a 
p x p matrix acting on the p states of the system. The 
different matrices corresponding to different possible ex- 
changes need not commute; hence, this type of statistics 
is called non-Abelian (braiding) statistics. Chern-Simons 
gauge theories with non-Abelian groups generically give 
rise to particles with such braiding properties. The lead- 
ing candidate to describe the v — 5/2 fractional quan- 
tum Hall stated is a state with non-Abelian (braiding) 
statistics2*i2*ii*i2ii 3 .. Other quantum Hall states observed 
in the first excited Landau leve&i^ might also be non- 
Abelian, possibly described by some of the states pro- 
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posed in refs. Il5lll6lll7l 

The different n particle trajectories which begin and 
end at the same positions, up to exchanges, are classi- 
fied by the elements of the braid group. If we do not al- 
low exchanges, as in the case of distinguishable particles, 
then we have the 'pure' braid group. Different varieties 
of particle braiding statistics correspond to different rep- 
resentations of the braid group. These representations 
are realized in the Hilbert spaces of many-particle sys- 
tems which support excitations with non-trivial braid- 
ing statistics and also in the corresponding Chern-Simons 
field theories. 

The Chern-Simons field theories which describe exotic 
braiding statistics are topological quantum field theories 
(TQFTs). The gauge fields in these theories do not probe 
local geometry - in particular, they disregard distance - 
so they only respond to topological properties, of the 
particle trajectories and also of the manifold on which 
they play out. Consider the amplitude associated with 
a process in which two pairs of statistics e 7Tl ^ m anyons 
and their anti-particles are created out of the vacuum. 
If the anyons braid around each other before the two 
pairs annihilate, the amplitude acquires a phase e 27rl£ / m , 
where I is the linking number of the trajectories. The 
non-Abelian Chern-Simons theories are related to more 
interesting link invariants such as the Jones polynomial. 
Even though the topology of the physical realizations 
which are probed in experiments is usually trivial, it 
is useful to consider more complicated topologies as a 
gedanken experiment probe of the underlying structure 
of the quasiparticle braiding properties. For instance, 
a system of statistics e 11 '™ anyons has a ground state 
degeneracy of m 9 on a surface of genus g. The phys- 
ical interpretation of these states is that quasiparticles 
pick up a phase e 27rm / m , n = 0, . . . ,m — 1 upon encir- 
cling the longitude of the torus. The underlying electrons 
have periodic boundary conditions in all of these states. 
By taking linear combinations of these states (or, in the 
language of conformal field theory, by applying the mod- 
ular S'-matrix), we can construct the states in which the 
quasiparticles acquire these phases upon encircling the 
meridian. 

Chern-Simons theories give a local description, but a 
redundant one. The basic idea is familiar in the con- 
text of electromagnetism: one could eschew the gauge 
field = (ip, A) in favor of the physically-measurable 
gauge-invariant electric and magnetic fields E, B but only 
at the cost of introducing non-local interactions between 
fields and charged particles in order to implement the 
Aharonov-Bohm effect. 

Indeed, any gauge theory will give rise to non-trivial 
braiding, by the Aharonov-Bohm effect. However, gauge 
theories with a continuous gauge group which are gov- 
erned by a Maxwell term alone are strongly fluctuating. 
In 2 + 1 dimensions, they fluctuate so strongly that they 
always confine charged (under the gauge group) parti- 
cles. Thus, the fractional excitations, which would have 
exhibited non-trivial braiding statistics if set free, are 



not, in the final analysis, part of the particle spectrum 
of the theory. A Chern-Simons term precludes such wild 
fluctuations by enslaving fluxes to charges. The other ex- 
ception occurs when the gauge group is discrete, as can 
occur on a lattice or in the continuum when a continuous 
group is spontaneously broken to a discrete subgroup. 
Discrete gauge fields have a deconfined phase, in which 
charges and fluxes (which are restricted to a discrete al- 
lowed spectrum) interact through the Aharonov-Bohm 
effect. 

This elegant field-theoretic description of exotic braid- 
ing statistics begs the question: are any of these possi- 
bilities actually realized in the real world? The answer 
is almost surely in the affirmative, at least for the case 
of anyons: the quasiparticle and quasihole excitations of 
fractional quantum Hall ground states are, according to 
theory, anyons. However, there has not, to date, been 
a direct observation of the phases associated with braid- 
ing. Beyond this, relatively little is known. In order 
to bridge this gap, it is important to understand in what 
other contexts one can observe topological phasesi^^ (or, 
equivalently, fractionalized phases) described by TQFTs. 
In particular, it is still an open question what types of 
microscopic Hamiltonians have, in their infrared limits, 
topological phases. 

We can derive some insight into this question by under- 
standing the structure of the Hilbert spaces of TQFTs. 
The Hilbert space of a TQFT in a topologically triv- 
ial geometry with no quasiparticles present is completely 
trivial: there is just a single state in the theory. However, 
when quasiparticles are introduced or the theory is put 
on a higher-genus surface, the Hilbert space will consist 
of a set of degenerate ground states. (In the former case, 
these are ground states for a fixed quasiparticle num- 
ber, assuming that the quasiparticles are not allowed to 
move.) The functional integral formulation of the theory 
encompasses all of these possibilities: one must simply 
perform the functional integral over fields defined on dif- 
ferent manifolds and with different boundary conditions. 
From a canonical point of view, however, the different 
possible topologies seem to give rise to completely differ- 
ent Hilbert spaces; they aren't even the same size, unlike 
in most systems, where changing the boundary condi- 
tions has little effect in the thermodynamic limit. Thus, 
it would seem that the canonical formalism is ill-suited 
for TQFTs. 

This is not the case. In fact, all of the Hilbert spaces 
associated with different topologies can be brought un- 
der the aegis of a mathematical structure called a mod- 
ular functot^. The modular functor can be enhanced 
by including edge states and an 'annulus category' which 
transforms these states. As we will discuss in this pa- 
per, this complex can be understood in more conven- 
tional physics terms as the representation theory of the 
commutator algebra of the fundamental gauge-invariant 
variables of the theory. The structure which is thus re- 
vealed gives, we believe, important clues about the nec- 
essary structure in any microscopic model which could 
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give rise to a topological phase. 

An important step in this direction is taken by the 
combinatorial construction of TQFTs. In this construc- 
tion, one builds the Hilbert space of a TQFT in the fol- 
lowing abstract fashion. Rather than start with the space 
of gauge fields A M , or some gauge-fixed version thereof, 
and complex-valued wavefunctionals on this configura- 
tion space, one begins with a set of 'pictures', collections 
of non-intersecting curves on a given surface (more pre- 
cisely, 1-manifolds). Two pictures are considered to be 
equivalent if they can be continuously deformed into each 
other. One then looks at the space of complex-valued 
functions on this set; it is an infinite-dimensional vector 
space. However, one can imagine imposing constraints on 
this vector space to reduce it to a finite-dimensional one. 
It turns out there is a limited number of ways of doing 
this, in fact just a single infinite sequence. If the wrong 
constraints are chosen, the resulting vector space will be 
zero-dimensional. Those favored conditions which lead to 
finite-dimensional vector spaces can be solved combina- 
torially. The construction can be generalized to surfaces 
with boundaries and 'punctures', at which curves can ter- 
minate. The latter are quasiparticles, and their statistics 
can be calculated by taking one puncture around another 
and using the constraints to simplify the resulting pic- 
ture. Such a representation of the states of a system 
in terms of loops will be familiar to some readers from 
analyses of dimer models, whose 'transition graphs' are 
configurations of loops^i. The TQFTs which we discuss 
are natural generalizations of the Z 2 gauge theory which 
emerges in the quantum dimer model on non-bipartite 
lattices^ 3 .. 

From the dimensions of their Hilbert spaces, we can 
guess that these TQFTs are closely related to but are not 
quite the same as known Chern-Simons theories. Unlike 
the latter, they are P, T- invariant. This comes about, in 
most cases, in an almost trivial way: the combinatorial 
construction leads to two decoupled Chern-Simons the- 
ories which are identical except that their chiralities are 
opposite. This theory is called the 'doubled' theory. The 
Hilbert space of the doubled theory is the tensor prod- 
uct of the Hilbert space of one copy of the Chern-Simons 
theory with its conjugate. As a result, the dimensions 
of the Hilbert spaces are the squares of the dimensions 
of the P, T-violating Chern-Simons theories. (In some 
'pathological' cases, the doubled theory is not simply the 
tensor square of the chiral theory but actually includes 
extra structure automatically repairing 'flaws' in the chi- 
ral theory.) The combinatorial approach automatically 
leads to the doubled theory: no artificial doubling of the 
degrees of freedom of the system was introduced at the 
outset; doubling emerged in the topological phase. 

These Chern-Simons theories are associated with 
gauge group SU(2). This is a consequence of the wonder- 
ful identification between SU(2) representations and di- 
agrams of non-intersecting curves (unoriented embedded 
1-manifolds) in a diskSi established in the Rumer-Teller- 
Weyl theorem. In short, SU(2) emerges from something 



very commonplace: sets of loops, which could represent 
domain walls, dimers, etc.. Z 2 or U(l) gauge groups 
occur in special cases in which they happen to coincide 
with a level k = 1 SU(2) theory, as we discuss further 
in section IV"Cl (The SU{2) gauge group is completely 
unrelated to the SU(2) spin-rotational symmetry which 
may or may not be preserved by a correlated electron or 
spin system which is in such a topological phase.) 

The combinatorial construction seems so far-removed 
from the ordinary methods of quantum field theory and 
many-body physics that it is initially somewhat surpris- 
ing that it leads to familiar TQFTs. However, a con- 
nection can be made by viewing these Hilbert spaces 
as representations of the commutator algebra of gauge- 
invariant Wilson loop operators rather than the less tan- 
gible gauge fields themselves. A simple representation of 
the form described above is only available for the doubled 
theory. The elucidation of this path to the combinatorial 
construction of the doubled theories is one of the main 
goals of this paper. 

The wavefunctions of these topological field theories 
are related to correlation functions in two-dimensional 
field theories, as is familiar from Laughlin's plasma 
analogy^ for quantum Hall wavefunctions and from the 
relation between Chern-Simons theory and conformal 
field theory 9 - 2 ^. 

In this paper, we will explore the physics of topo- 
logical phases observed with successively finer micro- 
scopes. First, we will discuss the long-wavelength ef- 
fective field theories which encapsulate the phenomena 
of exotic braiding statistics: Chern-Simons gauge theo- 
ries and discrete gauge theories. The former have ap- 
peared in the context of both the quantum Hall effect 
and also anyon superconductivity 28,29 The latter have re- 
cently been studied^S* 3 ^ in a revival of earlier ideas about 
RVB states 22 and spin-charge separation as a mechanism 
for superconductivit3r 21 i 32 i 33 i 34 i 35 i 36 . These two theories 
are the first in a sequence which forms the subject of this 
paper. The other members of this sequence are our main 
interest because they have excitations exhibiting non- 
Abelian braiding statistics. One of the principal motiva- 
tions for this work is the construction of models for topo- 
logical quantum computatiori 3 ^ 3 ^ 3 ^^^^; non-Abelian 
statistics (in fact, sufficiently rich non-Abelian statis- 
tics) is necessary in order to effect universal computa- 
tion purely through braiding operations because Abelian 
statistics results in the mere accrual of phases. Our ini- 
tial description is at the level of effective field theories, 
where one is at scales much larger than the characteris- 
tic size of quasiparticles, which are treated as point-like. 
By canonically quantizing these effective field theories, 
we can rephrase them in combinatorial terms. The com- 
binatorial construction of TQFTs has a natural inter- 
pretation as an intermediate scale description. At this 
scale, the system can be described in terms of fluctuat- 
ing curves. These curves must arise from shorter-distance 
physics and their dynamics must effectively impose cer- 
tain constraints which lead to the desired effective field 
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theories at longer wavelengths. The final step in this 
program is constructing microscopic models of interact- 
ing electrons, which will be the subject of a later paper. 
There, we will explore models which describe physics at 
the lattice scale. In order to produce the physics in which 
we are interested, these models must give rise, through 
their local interactions, to domain- walls or similar struc- 
tures which are the input for the combinatorial construc- 
tion. 

In essence, the current state of understanding is that 
topological phases are related to systems which, at low- 
energies can be understood as composed of fluctuating 
loops. In mathematical terms, the modular functor is lo- 
calized to curves - i.e. can be defined in terms of local 
rules on curves. This is, for the most part, the phe- 
nomenon which we describe in this paper. The next step 
is to find local Hamiltonians from which fluctuating loops 
emerge as low-energy degrees of freedom. This would be 
the localization of these modular functors to points. 



II. 'DOUBLED' CHERN-SIMONS THEORIES 

Chern-Simons theories break parity and time-reversal 
symmetries: a clockwise braid is not the same as a 
counter-clockwise one; it is its conjugate. There is a sim- 
ple way to start with such a theory and make a parity 
and time-reversal invariant theory: take two decoupled 
copies with opposite chirality. The resulting theory is 
called the 'doubled' theory or the 'squared' theory. The 
particle spectrum is now squared in size since it is the ten- 
sor product of the spectra of the two theories (which are 
mirror images of each other). This is a completely trivial 
procedure for a Chern-Simons theory, which hardly war- 
rants a special name. However, there are certain patho- 
logical topological field theories which can be fixed by 
'doubling' them and then specifying non-trivial braiding 
statistics between particles in the two tensor factors. As 
Drinfeld showed, this can be done in such a way as to 
cure the pathologies^. 

The Chern-Simons theories of the Halperin (m, to, to) 
quantum Hall states^ are examples of such 'pathologi- 
cal' theories. These are U(l) x U(l) Chern-Simons the- 
ories which can be written in a form in which one of 
the Chern-Simons terms has vanishing coefficient. Neu- 
tral excitations with non-vanishing spin quantum num- 
bers (or isospin if the second component is a second layer 
rather than the opposite spin) have trivial braiding prop- 
erties, unlike in the (to, m, n), to 7^ n states. As a result, 
the spin sector can condense and be gapless^, which is a 
'pathology' from the point of view of topological field the- 
ory. (From a more conventional physical standpoint, it's 
not a bug, it's a feature - and a remarkable one. This 
gapless excitation is a Goldstone mode associated with 
spontaneous breaking of the U(l) spin symmetry of an 
XY magnet.) It can be 'fixed' by Drinfeld's procedure if 
one introduces a second, oppositely directed, (to, to, to) 
state (which one might wish to call a (— m, —to, — to) 



state) and assigning relative statistics between the two 
spin sectors. It is not obvious that this can be done con- 
sistently in all such 'pathological' cases, but it can, as 
Drinfeld showed. 



'Doubling' initially appears very unnatural from a 
physical point of view. Why should two opposite chiral- 
ity copies of a theory arise? One might imagine sermons 
living on one sublattice of a bipartite lattice and anti- 
semions living on the other or, perhaps, up-spin semions 
and down-spin anti-semions. However, these sound more 
like clutching at straws in an effort to preserve time- 
reversal invariance; in cither case, the combination of 
time-reversal with translation or spin-rotation is still bro- 
ken. However, 'doubling' can occur in a more organic 
way. In fact, it is not just a completely natural occurrence 
in our models, it is absolutely unavoidable. For instance, 
the simple model given in section lVIIIl is a 'doubled chiral 
spin-liquid' while, thus far, no model with short-ranged 
interactions has been shown to have an undoubled chi- 
ral spin-liquid ground state. In order to eliminate dou- 
bling, additional structure - possibly unnatural - must 
be added to the models. 



From an algebraic point of view, it is useful to observe 
that in a system in a magnetic field, translations in the 
x and y directions don't commute. This is expressed in 
Chern-Simons theory by the statement that in the canon- 
ical formalism, the x and y components of the gauge field 
do not commute. Thus, they cannot be simultaneously 
diagonalized. In this sense, a Chern-Simons gauge field is 
'half of a gauge field: only half of its degrees of freedom 
can be specified in a given state. By doubling the theory, 
or taking two copies, we have, in a sense, doubled 'half 
of a gauge field, thereby yielding one gauge field. Thus, 
it is possible, in the doubled theory to have a basis of 
states which are eigenstates of the Wilson loop operators 
- which are the natural gauge-invariant operators - as- 
sociated with a gauge field. In the Abelian case, this can 
be done completely explicitly. We can construct a gauge 
field from the two opposite chirality gauge fields, and in 
its diagonal basis, the Wilson loop operators are analo- 
gous to number operators N, and operators which do not 
commute with them, such as the Wilson loop operators 
for either of the chiral gauge fields, are roughly analo- 
gous to creation and annihilation operators , a. In the 
non- Abelian case, this basis cannot be related in such a 
simple way to the original gauge fields, but roughly the 
same structure is present. 



We will see that systems which have configurations 
which can be mapped onto those of a loop gas have a 
Hilbert space which naturally admits such a basis. Thus, 
if they enter a topological phase, we expect it to auto- 
matically be a doubled one. 
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III. PARTON MODEL CONSTRUCTION FOR 
CORRELATED ELECTRON SYSTEMS 

A useful heuristic for understanding the physics of 
a correlated electron Hamiltonian (in zero or non-zero 
magnetic field) involves rewriting the electron opera- 
tor in terms of auxiliary 'parton' operators or slave 
fermion/boson operators. This strategy can also be ap- 
plied to interacting bosons models, which could be real- 
ized in materials in which Cooper pairing leads to the 
emergence of effective bosonic degrees of freedom corre- 
sponding to the pairs. In this section, we will use the 
parton formalism to describe the phases of interest. 

Consider a model of spins on a lattice. Suppose that 
there is spin S = N/2 at each site. The Hamiltonian may 
be of the form 



H = J £ S t ■ S j 



(1) 



<i,3> 



where the ellipses denotes other terms including, perhaps, 
ring exchange terms or next-neighbor interactions. We 
now introduce an SU(N) multiplet of spin-1/2 fcrmions 
faa, with a = 1,2, ... ,N and a =|, J, so that 



Pi 



(2) 



The fermions must satisfy the constraints that there be 
N fermions per site: 



faaifaai — N 

and that there be a color singlet at each site 

faai^abfapi = 



(3) 



(4) 



k 9 

where T~ b , k = 1,2,..., N — 1 are the generators of 
SU(N) in the fundamental representation. These con- 
straints guarantee that there is spin N/2 at each site. 

Using this representation, we can rewrite the Hamilto- 
nian in the form 



H — J faaifbajfbpjfa. 



0i 



(5) 



<hj> 



In mean-field approximation, this Hamiltonian can be 
written as: 



h = -j E u *b,vfLfi 



baj 



(6) 



<hj> 



where 



U ab - (fbpjfaPi) ( 7 ) 

Consider the following mean-field solution: 



Ull = (jf) e*« S ab 



(8) 



where i c ff is a parameter to be determined from the 
saddle-point condition. Assuming that the system is on 
a square lattice of side a, let us suppose that 



77 i 



x i) (Vj +Vi) I * 



so that 



n 

plaq. 



= e 2TTi/k 



(9) 



(10) 



Here, k is assumed to be an even integer. In this solution, 
the system generates an effective magnetic field for the 
/ aa s such that there is half of a flux quantum per pla- 
quette. The commensurability of the field ensures that 
precisely k/2 Landau levels are filled by each of the / oa s. 
(Landau levels are broadened into Hofstadter bands since 
we are on a lattice. The fermions fill an integer number 
of these bands.) 

Does this solution actually occur? Is it stable, i.e. is 
it a local minimum of the action for this system? The 
answers to these questions depend, for the most part, 
on the particulars of the given Hamiltonian. Later in 
this paper, we will take some steps towards addressing 
them. While we can't determine the absolute stabil- 
ity of this solutions using the approach of this section, 
we can test its local stability against small fluctuations. 
There are certain fluctuations which, a priori, are likely 
to be important. Consider the gauge transformations 
/«« -> KbhauUZ - W^U&Wg, where W* b are 
U(N) matrices assigned to the points i of the lattice. 
These transformations leave the Hamiltonian invariant. 
Hence, the class of mean-field configurations 



IP3 = I ^ 
aC [ t 



J<Pi. 



(11) 



have vanishing energy cost if gij = Aj — Aj for some u(N) 
Lie algebra- valued matrices Aj, i.e. if is pure gauge. 
The energy cost should be a gauge-invariant function of 
gij] since gauge field fluctuations can be large we must 
consider them carefully. Let us break the U(N) gauge 
field g^ into U(l) and SU(N) parts, Cy and a^, respec- 
tively. The latter is a traceless, Hermitian, NxN matrix, 
i.e. an SU (N) Lie algebra- valued field. The time compo- 
nents of these gauge fields, Co and ao, can be introduced 
as Lagrange multipliers which enforce the constraints © 
and 10} . 

To see that the gauge field fluctuations do not desta- 
bilize our mean-field solution, we integrate out the / aa s. 
The crucial point is that the f aa s fill k/2 Landau levels 
in the mean-field solution. By standard arguments, we 
then have the following effective action for c and a 



S[a] = NkS cs [c} + kScs[a] 
Nk f 
= — I e'" •-<■„</„<■„ 
An , 
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Nk 

47T 



c A dc 



k f ( ,2 
-— / tr a A da + -a A a A a 
Air I V 3 



(12) 



The effective actions for the U(l) gauge field c M and the 
SU{N) gauge field af = a% T° are the corresponding 
Chern-Simons actions at level k, i.e. with coupling con- 
stant k. fijk are the structure constants of SU (N). Since 
each f aa fills k/2 Landau levels, their response to an ex- 
ternal gauge field must break P, T and be proportional 
to k; the requirement of gauge invariance then dictates 
(|12J) . up to multiplication by an arbitrary integer. A di- 
rect calculation shows that this integer is 1. For the U(l) 
gauge field, it is just the Hall conductivity of k/2 filled 
Landau levels of both spins and N colors. We have sup- 
pressed subleading terms, such as the Maxwell terms for 



tain types of models, but, for now, let us simply consider 
'doubled' theories as another logical possibility. 

Let us consider a model in which there is an integer 
spin N at each site of the lattice. Suppose we now intro- 
duce 2N spin— 1/2 fermions f a Aa, with a = 1, 2, . . . , JV, 
A = ±, and a =|, J, so that 



Si — flAai a apfaAl3i 



(13) 



The fermions must satisfy the constraints that there be 
2N fermions per site: 

fLfaai = 2iV (14) 

and that there be a color singlet at each site 

flAai T abABfaBi3i = (15) 



and af, which should appear in eq. d since they are where T~ bAB , k = 1,2,..., AN 2 - 1 are the generators 



irrelevant in the low-energy limit. 

The evenness of the level k results from the presence 
of equal densities of up- and down-spin fermions. If spin- 
rotational symmetry were broken by the presence of dif- 
ferent densities of f a A]S and / a A|S, then odd level k could 
also result. 

The effective field theory ljl2|l is a gapped theory, so 
gauge field fluctuations about this mean-field solution are 
not large 45 . In other words, if the configuration (JSJ is 
a saddle-point for some Hamiltonian, then this saddle- 
point is stable against fluctuations of the form ljll|l . Said 
differently, the breaking of P, T permits the appearance 
of a 'mass term' for the gauge field - the Chern-Simons 
term - which renders the phase stable. 

The appearance of a Chern-Simons term has another 
important effect: the excitations of the theory have exotic 
statistics. As a result of their interaction with the gauge 
field, the f aa s are anyons and, (except for the special 
case k = 1 which cannot occur anyway in this construc- 
tion), their braiding statistics is non-Abelian. (It is not 
quite consistent to talk about the / aa s interacting with 
a Chern-Simons gauge field ajt since we had to integrate 
out the f aa s in order to generate the Chern-Simons term 
for aj,,. However, we can introduce external source fields 
for the f aa s. When the / aQ s are integrated out, a cou- 
pling between the source fields and ajt is generated, so 
that when these sources are braided, the result is non- 
trivial.) 

The phases which we have just discussed share the 
attractive features that they are stable and that they 
support quasiparticle excitations with exotic braiding 
statistics. They also spontaneously break parity and 
time-reversal invariance. The latter would seem, at first 
glance, to be a necessary condition for exotic statistics. 
However, the closely-related 'doubled' theories, which we 
discussed in the previous section, also support exotic 
statistics, but they do not break P, T. This is a useful 
feature in, for instance, a theory of a material in which 
P, T violation has been experimentally ruled out. We will 
later discuss why such phases are likely to occur in cer- 



bABi 

of SU(2N) in the fundamental representation. These 
constraints guarantee that there is spin N at each site. 

Following the steps which we made earlier, we make a 
mean-field approximation in which the Hamiltonian takes 
the form: 



H 



J U a bAB,ijf a A ai fbBaj 
<i,3> 



where 



u abAB 



, fbBPjfaA/3i 

and consider the following mean-field solution: 



(16) 



(17) 



U ab++ 




) S, 


Tjij 
U ab 




j e -i4>i) 






Tjij 
U ab+- 


- u ij 

" U ab- + 


= 



with 



(18) 



(19) 



Again, k is assumed to be an even integer. The system 
generates an effective magnetic field for the / act s in which 
k/2 Landau levels are filled by each of the f aa s. 

Turning now to fluctuations about this saddle-point, 
we see that there are two SU (N) gauge fields to go with 
the U(l) gauge field. The low-energy fluctuations are of 
the form 



U ab^ 



ab — 
u ab+- 



tcS 

T 

tcft 

J 



e v 3 e ■ 



-i<t>i) picij 



(e*4) 

V J ab 
\ J ab 



— u lJ — n 

— U ab — h ~ U 



(20) 



The theory actually has a U(\) x SU(2N) gauge sym- 
metry. However, the saddle-point is only invariant under 
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global 17(1) x SU(N) x SU(N) transformations. The 
other global SU(2N) transformations are broken at the 
saddle-point level and, hence, the corresponding gauge 
fields are massive by the Anderson- Higgs mechanism. 

Equivalently 35 , symmetry allows a term in the effective 
action of the form 



tr (Fi 



ijkl ^ 



'Fi 



i' j'k'V & 



(21) 



where ijkl and i' j'k'V are two plaquettes which are one 
lattice spacing apart so that they are connected by the 
link ji' . Fijki is the product of the [7's around the pla- 
quette ijkl. At the saddle-point, Fiju = fVj'fc'i' = F = 
diag(e 27 ™/ fc , e~ 27 "/ fc ) which commutes with the elements 
of SU(N) x SU(N) but not the remaining elements of 
SU(2N). Expanding (|21[l to second order in ay^ we find 
a term 



tr {F\[F,a iV \,a iV \) 



(22) 



which is a mass term for the gauge fields associated 
with SU(2N) generators which don't commute with 
diag(e 27ri / fc ,e- 27 "/ fe ), i.e. those not in SU(N) x SU{N). 

The effective action for the U(l) x SU(N) x SU{N) 
gauge fields can be derived by integrating out the 
fermions. We find: 



S[a] = k S C s [a + ] - k Scs [a 



(23) 



Since the f a — a fermions move in an effective field of flux 
— 2n/k, there is a negative sign in front of the second term 
in (I23JI . Furthermore, there is a cancellation between 
the contributions of the f a + a and f a - a fermions to the 
coefficient of the Chern-Simons term for the U(l) gauge 
field c M . Hence, we have two opposite chirality SU(N) 
Chern-Simons gauge fields. They are gapped and lead 
to exotic braiding statistics, but they preserve P and T 
since we can exchange a + and a~ when we perform time- 
reversal or a parity transformation. We will call these 
theories SU(2)k x SU(2)k Chern-Simons theory, with the 
overline signifying that the second Chern-Simons term 
has opposite chirality. 

The U(l) gauge field does not have a Chern-Simons 
term. Thus, it is in one of two massive phases. (1) A 
confining phase in which gauge field fluctuations cause 
[7(1) charge to be confined and c M to be gapped. (2) A Z 2 
phase which results when eab^AB^apfaAafbBp condenses, 
breaking C7(l) to Z 2 . The remaining Z 2 gauge field has 
a phase in which Z 2 charges are deconfined. The former 
theory can be considered as a subset of the latter: both 
of these phases have quasiparticles with exotic braiding 
statistics, but the former only has those which are created 
by fermion bilinears such as f aAa fbBp while the latter has 
the full set of fractionalized quasiparticles. In the SU (2) 
case, the former theory has only integer-'spin' (i.e. under 
the internal SU(2) x SU{2) gauge symmetry) quasiparti- 
cles while the latter has half-integer as w ell. We will call 
the former the 'even part' of SU(2)k x SU(2) k Chern- 
Simons theory. 



We have introduced these phases in the context of spin 
systems, but they can arise in a variety of contexts. The 
above analysis can be extended to finite doping in the 
usual way. We can also expect such phases in, for in- 
stance, correlated boson models. Consider an example of 
the latter, 



H = —t 



<i,j> 



J2 VijNiNj (24) 



<i,3> 



where i,j label sites on a lattice, B>i creates a Cooper pair 
on site i, and Ni = -Bj-Bj. We will assume that Vu = oo 
so that we have hard-core bosons. Let us assume that 
there are precisely half as many bosons as lattice sites. 
We assume that Vij > 0. For Vij small, the system will 
be superconducting. For Vij sufficiently large, the sys- 
tem will be insulating. We would like to explore possible 
insulating phases. In order to do this, we write B>i in the 
following representation: 



Bi = 



1 



fNl 



fa\i fa 



Jap 



(25) 



where = 1, 2, . . . , N . We have introduced an even 
number N of auxiliary fermions j a i such that the boson 
Bi is an SU (N) symmetric bound state of them. 

Focusing on the hopping term, we see that the Hamil- 
tonian can now be written as: 

H = ~ ~^7j ^""^ { e a,ia. 2 ---a N ft 1 ifa 2 i---fi N i X 
' <i,j> ^ 

e b t b 2 ...b N fbxj fb 2 j ■ ■ ■ fb N j + h.C.^ + . . . (26) 

The ellipses represent next-nearest-neighbor and longer- 
ranged interactions. In mean-field approximation, this 
Hamiltonian can be written as: 



(27) 



<hj> 



where 



7^aa 2 ...ajv^bb 2 ...b]v ^ 



^/aii ■ ' ■ fa N i fb 2 j ■ ■ ■ /6jviy(28) 

Consider the following mean-field solution: 

tfj b ={^f)e i ^8 ab (29) 

where 

<t>ij = -^f (xj - Xi) (yj + yi) /a 2 (30) 

for some integer m. Only such fluxes are allowed since 
the bosons B must see vanishing flux. If the / a s see 
flux 2-KmjN then the Bs see flux 2itm which is gauge 
equivalent to zero. Suppose that 2m is a divisor of ./V 
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so that k — N/2m is an integer. Then the / a s fill k 
Landau levels and, integrating out the fermions, we find 
that the effective field theory for fluctuations about this 
saddle-point is SU(2mk)k Chern-Simons theory. 

The doubled version of this theory can be obtained by 
using a representation with two sets of ./V fermions f +a 
and f- a : 



with a mean-field 




l a 2 ...a N f+aii f+a 2 i ■ ■ ■ f+a N i 
e aia2...ajy f—aii f—a,2i • • ' f—aNi (^-0 



IT 13 

U ab+ + 



U ab- + 







(32) 



The effective fie ld theory fo r fluctuations about this state 
is SU(2mk)k x SU{2mk) k . 

This construction can be adapted to a system of itin- 
erant electrons, for which we introduce an odd number 
of partons whose bound state is an electron. In this way 
and also by using other straightforward generalizations 
of the above constructions, we can obtain other phases 
of correlat ed electro n systems which are described by the 
SU{N) k x SU(N) k family of effective field theories. 

In the remainder of this p aper, w e will focus on the 
case N — 2, i.e. the SU(2) k x SU(2) k Chern-Simons the- 
ories. These can actually arise in two different ways from 
the above constructions. One is via the obvious path: a 
spin-2 magnet, for instance, can have an SU (2) k x SU(2) k 
phase. The alternative route follows from rank-leve l dual- 
ity: a spin-fc magnet can have an SU(k)2 x SU(k)2 phase, 
which is equivalent to SU(2) k x SU(2) k . The latter con- 
struction suggests that the higher-level phases should be 
sought in higher-spin magnets, an observation on which 
we will comment further at the end of this paper. 

The parton or slave-particle representations which we 
have discussed in this section have a Z 2 gauge symmetry 
under which /+ and /_ are exchanged independently at 
each point. (In the spin systems which we discussed ear- 
lier, this Z 2 is a subgroup of SU(2N).) This symmetry is 
broken at the saddle-point. P, T are also broken. How- 
ever, the combination of a Zi gauge transformation and 
T or P leaves the saddle-point invariant. Thus, we have 
the symmetry- breaking pattern Z?^ uge x — > Z 2 iasona . 
In this way, P, T are preserved since they can be identi- 
fied with the diagonal Zi which remains. This is similar 
to the preservation of rotational invariance in monopole 
solutions of gauge theories: the solutions are not rota- 
tionally invariant, but the combination of a rotation and 
a gauge transformation leaves the solution invariant. 

To summarize, we have seen in this section how the 
parton or slave-particle formalisms often used to discuss 
strongly-correlated electron systems can lead to stable 



phases corresponding to doubled Chern-Simons theories. 
Thus far, we have not discussed why the doubled theo- 
ries might arise rather than the undoubled ones or even 
some entirely different phases. In other words, the ques- 
tions which we have not addressed is why should one 
of these mean-field solutions have the minimum energy? 
In order to begin to answer this question, we will take 
a somewhat circuitous route which will involve solving 
them first. Special properties of the doubled theories 
will emerge once we have discussed them more fully. In 
order to understand the physics of these theories, we will 
need a set of non-perturbative methods, discussed in the 
following sections. 



IV. SOLUTION OF DOUBLED ABELIAN 
THEORIES 



The Chern-Simons theories (both undoubled and dou- 
bled) which we have just encountered are topological field 
theories. Thus, they have only a finite number of degrees 
of freedom and are 'trivial' in the sense that they are 
a problem in quantum mechanics, rather than quantum 
field theory. On the other hand, the physics which they 
describe is non-trivial precisely because it is topologically 
invariant and measurable physical properties probe the 
topology of the manifold on which they live. Thus, these 
theories are soluble, but not with the methods ordinarily 
used to solve field theories. In particular, we will need 
to use the canonical formalism almost entirely, although 
ideas from the functional integral approach will prove 
useful. In this section, we present the solution of these 
theories, which takes the form of a construction of their 
Hilbert spaces, together with physical observables acting 
on these Hilbert spaces. Their algebraic structure re- 
veals connections with theories of exactly-soluble statis- 
tical mechanical models and, perhaps more importantly, 
offers some valuable clues about which microscopic mod- 
els might give rise to these topological phases. In this 
paper, we will make a few comments about microscopic 
models and defer a more serious discussion to a following 
paper. 

Topological field theories have vanishing Hamiltonian, 
so they do not describe dynamics in the ordinary sense in 
which, say, a non-linear sigma model would. They only 
describe the braiding properties of quasiparticles and the 
sensitivity of the system to the topology of the surface 
on which it is realized. Consequently, the entire problem 
of solving these theories amounts to constructing their 
Hilbert spaces; there is no energy spectrum to compute 
because all states have zero energy. In other words, the 
problem is an algebraic one of constructing these Hilbert 
spaces as the representation spaces of the commutator 
algebra of the physical observables of the theories. In 
this section, we solve this problem for the 'warm-up' case 
of Abclian theories. 
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A. Abelian Chern-Simons Theory 

Let us begin our discussion of the effective field theo- 
ries for topological phases of interacting electron systems 
with one of the simplest such theories, Abelian Chern- 
Simons theory: 



5 : 



in 
Air 



(33) 



The coefficient is chosen to be an integer divided by Air in 
order to ensure invariance under large gauge transforma- 
tions. A more refined effective field theory could include 
further terms such as a Maxwell term. However, these 
terms are irrelevant in the infrared, so we drop them. 

Let us assume that we are working at energies much 
lower than that required to excite any of the particles to 
which the gauge field is coupled, so that we can consider 
(|33|l in isolation. This theory would appear to be com- 
pletely trivial. Suppose we take Coulomb gauge, a = 0. 
In taking this gauge, we must remember to impose the 
constraint which follows from varying do in (|33|l 



V x a = 

In this gauge, the Lagrangian takes the form 

£ = ^-a 2 d Q ai = ITidoai - 

where the canonical momentum conjugate to a\ is 

dC m 



IIi 



<9(<9o fl i) 27r 



«2 



(34) 



(35) 



(36) 



From the second equality of 1)35(1 , we see that the Hamil- 
tonian vanishes. Thus, the effective action only describes 
the ground state - or states. The interesting structure 
of the theory follows from the canonical equal-time com- 
mutation relation 



[0l(x,t),n 1 (3C',*)] = jj(x-x') 



27TI 

[ai(x,i),a 2 (x',t)] = 5(x-x') 

m 



(37) 



(38) 



Note that the all-important integer m appears here in 
the commutation relation because it appears in front of 
the Chern-Simons action. 

On the infinite plane or the sphere, the ground state 
is a unique, non-degenerate state. Pure Chern-Simons 
theory (i.e. without any other fields coupled to it) has no 
other states. However, suppose that the theory is defined 
on the torus. Then a gauge field satisfying V x a = can 
still give rise to a non-trivial holonomy W[j] around the 
closed curve 7 if 7 winds around one of the non-trivial 
cycles of the torus. 



a i§ a-d\ 



According to the constraint, W[y] does not depend on 
the precise curve 7 but only on how many times it winds 
around the generators of the torus, i.e. on its homotopy 
class. Furthermore, it is clear that W[j] is multiplicative 
in the sense that its value for a curve 7 which winds twice 
around one of the generators of the torus is the square 
of its value for a curve 7 which winds once. Hence, we 
have only two independent variables. This is revealed by 
solving the constraint. In a coordinate system 81,62 € 
[0, 2tt] on the torus, we have a = (a±/2n, a 2 /2ir) + V<p. 
If we take ip to be single- valued, then invariance under 
large gauge transformations requires that we make the 
identification a,; = a,; + 1. Wi and a* are related by 



Wi 



(40) 



where gammai, i — 1,2 are representative curves which 
encircle, respectively, the meridian and longitude of the 
torus once. 

From (|38|) . we have the following equal-time commu- 
tation relations: 



r 1 - 2lT 

[Oil, Ct 2 \ = I 



(41) 



Since a±, a 2 are not themselves gauge- invariant, we can- 
not simply use the analogy between their commutation 
relations and those of p, x for a single particle. We must 
work with the gauge invariant quantities Wj. Since 



g-tCKi g ia 2 _ [qi,Q2]/2 e iui+iu2 



we have the commutation relation 



W X W 2 = e 27Ti/m W 2 W 1 



(42) 



(43) 



This algebra can be implemented on a vector space in 
the following way: 



Wi \n) 
W 2 \n) 



e 2Tmi/m | n ) 

ln + 1) 



(44) 



(39) 



with n G Z. This is a representation of this algebra on 
an infinite-dimensional Hilbert space. However, there are 
finite-dimensional truncations which also allow a repre- 
sentation of this algebra. Suppose that we simply restrict 
n G {ii 3 + 1> 3 + 2, • • • , j + qm — 1} for any integer j and 
any positive integer q. and define W 2 \j + qm — 1) = \j). 
Such a vector space has dimension qm and it gives a per- 
fectly good representation of the algebra. Which value of 
q gives the Hilbert space of U (l) m Chern-Simons theory? 

To answer this question, first note that the phase space 
of the classical theory has finite volume. It is parameter- 
ized by dip € [0, 2tt] defined in 1)41) [I. This is analogous 
to the case for a single spin, where phase space is the sur- 
face of a sphere S 2 , but it is unlike the case of a particle 
constrained to lie on a circle, where the coordinate takes 
values on a circle but the momentum is arbitrary, so that 
phase space is S 1 x R. From the commutation relation 
for ot\, a 2 , we see that ot\ and ma 2 /2ir are canonically 
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conjugate coordinates on phase space; thus phase space 
has volume 2irm. We expect the dimension of Hilbcrt 
space to be roughly equal to the volume of phase space, 
measured in units of h = 2ith = 2tt (since we have set 
% = 1). In the classical - or large m - limit, this should 
be an exact equality, so the only choice is the minimum 
possible one allowed by the algebra (|43|l . q = 1, i.e. di- 
mension to. 

We can restate this by saying that the Hilbert space of 
the theory is obtained from the infinite-dimensional one 
of H44fl by requiring that all physical states be annihilated 
by the projection operator 

J2(\n)-\n + m))((n\-(n + m\) (45) 

n 

To summarize, the Hilbert space of U(l) m Chern- 
Simons theory is spanned by the basis vectors \n) with 
n = 0, 1, ... ,?7i — 1, i.e. the ground state is m-fold degen- 
erate. On a genus g manifold, this generalizes to m 9 . The 
inner product {n\n') — § nn i is fixed by the requirement 
that Wi and W 2 be unitary. 

This has an interpretation in terms of the 
(quasi)particle spectrum of the theory - about which we 
might have thought that we would lose all information at 
low energies. Imagine creating a quasihole-quasiparticle 
pair, taking them around the meridian of the torus and 
annihilating them; call the corresponding operator T\. 
Let T2 be the operator for such a process around the 
longitude. If the quasiparticles have statistics ir/m, then 

TiT 2 = e 2 " /m T 2 Ti (46) 

because the particles wind around each other during such 
a process, as depicted on the right of figure This is 
precisely the same algebra <|43|) which we found above, 
with representations of minimal dimension m. 




FIG. 1: The operation T^T^T^ 1 results in a phase e 2 ™ /m 
because it is equivalent to the braiding operation on the right. 

Hence, if we know that the ground state degeneracy of 
a system on a genus-<? manifold is to 9 , then one expla- 
nation of this degeneracy is that it has non-trivial quasi- 
particles of statistics 0,7t/to, . . . , (to — l)7r/m. 

One awkward feature of the basis which we have just 
constructed for the Hilbert space of Abelian Chern- 
Simons theory is that we had to choose a particular di- 
rection on the torus. If Q\ is the coordinate along the 
meridian of the torus, then the Wilson loop operator W\ 
corresponding to the meridian is diagonal in this basis, 
while the Wilson loop operator taken along the longitude 
is not. This appears to be unavoidable because the two 
operators do not commute. However, by considering the 
'doubled' theory, a theory with two Chern-Simons fields 



with equal but opposite coupling constants, to and —to, 
we can have a more democratic-looking Hilbert space. 
This is more than just an aesthetic requirement, since 
the simplest microscopic models give rise to such Hilbert 
spaces, as we will see later. 

B. Doubled Abelian Chern-Simons Theory 

The Hilbert space of the theory with action 

s = 4^ J ^ Vp ^ v a p - — J e^P Cfl d u c p (47) 

is clearly just the tensor product of the Hilbert spaces 
associated with the two terms in the action or, in other 
words, the tensor product of the Hilbert space (|44|l with 
its complex conjugate (the conjugation results from the 
minus sign in front of the second term in eq. I47f) . We 
will usually call this theory U (l) m x U(l) Chern-Simons 
theory, but we will sometimes call the theory the 'double' 
or 'Drinfcld double' of U(l) m Chern-Simons theory. This 
terminology is unnecessary in this simple case, but there 
are instances of pathological chiral topological field the- 
ories whose pathologies can be cured by doubling them 
and adding some extra structure, as Drinfcld found. 

For simplicity, let us consider the case m = 2. Now, 
following l|39|) , we can define the operators W+ [7] associ- 
ated with the gauge field and the analogous operators 
W- [7] associated with c M . We now have the tensor prod- 
uct of two operator algebras: 

W + h]W + W] = (-l/^TT+M^+M 

W_[7]W_[Y] = (-l) /(7 ' 7 ° I4M 7 '] WM7] 
W+[j] wJ[i] = WJc/W+h] (48) 

I(a, 7) is the intersection number of a and 7. It is use- 
ful to define the operators L[y] = W+ [7] W- [7] , which 
commute with each other: 

L[ 1 ]L[ 1 '} = L[ 1 ']L{ 1 ] (49) 

We introduced these operators to emphasize the point 
that (|48|l is not the tensor product of just any two alge- 
bras pulled off the street, but of two identical ones. As a 
result, there is a commuting set of operators L[y] associ- 
ated with curves 7. For future reference, we display the 
relations which the L[y]s obey with W+fy], W_[7]: 

mW-\i] = (-l) /(7 ' y) WM 7 ']£[7] (50) 

One way of representing this algebra is on the Hilbert 
space |n+, 7i_) with 71+, n_ = 0, 1 

LiK,n_) - (-!)»+-«- 

L 2 |n+,n_) = I n+ + l,7i_ + 1) (51) 

where L1.2 = £[7] with 7 a meridian or longitude, re- 
spectively. W+ [7] and W- [7] act as one would anticipate 
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from (|44(l . As expected, the ground state degeneracy is 
m 2 = 4. 

Since L1.2 commute, they can be simultaneously diag- 
onalizcd. They are diagonal in the following basis 



space which serves as our pre-Hilbert space. 



\li,£i) = 52(-l) nta \n + £ u n) 



(52) 



n=0 



(The generalization to arbitrary m is clear: the sum 
ranges from to m — 1, and (—1) is replaced by e 27 ™/ m .) 
Then 



£i >2 |4,4> = (-ir' 2 |4,4) (53) 



Meanwhile 



W +1 \l u h) = (-1)^|4,4 + 1) 
W +2 \£ 1 ,£ 2 ) = |4 + 1,4) (54) 

with similar relations for W-. Thus, we can think of 
L1.2 as being analogous to number operators Nk while 
W+1.2 are analogous to raising operators a k . The anal- 
ogy is not quite right because W+i 2 do not commute 
with each other. However, this is a useful analogy never- 
theless. (On the torus, we actually can take a commuting 
set of 'raising/lowering' operators, W+i and W-2)- We 
have one such 'number' operator for each generator of 
the torus. 

Note that states |\&) can be written as wavefunctions 
in this basis: 



(55) 



*[4>4] maps two integers modulo 2 into the complex 
numbers. 

A slightly more abstract representation of the same 
Hilbert space will prove useful when we generalize this 
construction to more complicated theories. It will also 
suggest connections with microscopic models. Again, we 
will work on the torus, but the extension to other sur- 
faces is straightforward. The basic idea is to define wave- 
functions on the space of curve configurations on a given 
surface. The discussion of the previous paragraph can be 
framed in these terms if we think of t\ , I2 as defining a 
topological class of curves. 

In order to do this in more general terms, we need 
a few definitions. We define the following notation: let 
{a} be an isotopy class of one-dimensional submanifolds 
of the torus. A one-dimensional submanifold, a, of the 
torus is simply a set of non-intersecting smooth curves. If 
the submanifold is connected, then it is a curve; however, 
we want to allow the multi-component case. We will use 
the term multi-curve to denote such a one-dimensional 
manifold, a and ol are in the same isotopy class if they 
can be smoothly deformed into each other4& We define 
a 'pre-Hilbert space', H by associating an abstract vec- 
tor I {a}) to every isotopy class, {a}, of one-dimensional 
submanifolds of the torus. By forming all linear combi- 
nations with complex coefficients, we arrive at the vector 



Ha} 



e C 



(56) 



The vectors in our pre-Hilbert space are complex-valued 
functionals of the isotopy classes of one-dimensional sub- 
manifolds of the torusfe 



rp[{a}] = ({a}\i>) 



(57) 



We called this the 'pre-Hilbert space' of our theory 
because the actual Hilbert space of the theory, H., is a 
subspace of ft. We define Ti as the subspace consisting 
of ^[{ajjs satisfying the following constraints: 



*[{«}] = 
*[{"}] = 



-*[{aUO}] 
-¥[{&}] 



(58) 



aU O is the one-dimensional submanifold of the torus 
which is obtained from the union of a with a contractible 
loop; a is obtained from a by performing the cutting and 
rejoining operation )(— > ^ on any part of a (see fig. |5J. 
Since this operation, which is called 'surgery', can be 
performed on any part of a, the relation H58|) must hold 
for all possible {a}. 





FIG. 2: A surgery operation of this type changes the wave- 
function by a factor of —1. The shaded area represents an 
arbitrary multi-curve which completes the 1-manifold shown. 

Said differently, the Hilbert space of our theory is the 
subspace of pre-Hilbert space which is annihilated by the 
two projection operators: 

K-i = (lW) + |{aUO») ((Ml + ({a U Oil 

P2,-i = (iW) + !{«}>) ((WI + ({«}| X 



(59) 



As before, the first projector must be applied for all pos- 
sible contractible loops and the second projector must be 
applied for all possible {a}. 

This Hilbert space is four-dimensional on the torus. 
To see this, observe that the operations a — * a U O 
and )(— > <— - preserve modulo two the winding number 
of a multi-curve about the torus but do not preserve the 
winding number itself. There are four possible winding 
numbers about the meridian and longitude of the torus: 
(0, 0), (1, 0), (0, 1), (1, 1). We can define a function 
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which vanishes on all isotopy classes of one-dimensional 
submanifolds of the torus which do not have winding 
number (i,j) modulo two. We can assign ^a,j) the value 
one for some isotopy class which has winding number 
modulo two; it takes value ±1 on all other isotopy 
classes which have winding number (i,j) modulo two, 
according to the relation (|58(l . These four wavefunctions 
^(ij) form a basis of Hilbert space. A further point must 
be checked: that no 1-manifold a can be related back to 
itself by an odd number of surgeries. This would intro- 
duce the relation a = (— l) odd a or a — 0. This fact is 
intuitively obvious, but a proof would take us on a topo- 
logical digression, so we omit it here. 

Thus, our Hilbert space is four-dimensional, as it 
must be if it is to be the same as the Hilbert space 
of U(l)2 x U(l) 2 Chern-Simons theory, which we con- 
structed earlier. To see that it, indeed, fulfills its raison 
d'etre, namely to furnish a representation of the algebra 
(03 1 (EH) w e define the operators W + [7] , W" [7] : 

VK+[ 7 ] = i" (7 < Q) <F[{aU n7 }] 

- (-l)^> a > *[{«}] (60) 

rt(7, a) is the number of intersections between 7 and a 
(without regard to sign, unlike 1(7,0;)). The notation 
aU H 7 simply means the union of a and 7 if they do not 
intersect; if they do, then the crossing is resolved by turn- 
ing to the right as the intersection is approached along 
a, as shown in figure |3J aL) L j is defined in the oppo- 
site way. It may not be obvious that aU H 7 depends only 




FIG. 3: When a and 7 intersect, aUnj is defined as shown 
above. 

on the isotopy class of a. For instance, suppose that a 
and 7 do not intersect. Then, a can be continuously de- 
formed into isotopic a' such that a' intersects 7 twice. 
Thus, if we act on ^[{a}] with VF + [{7}], we would seem 
to get a different factor, j™(7>«) = l or i n ^- a ') = -1, 
depending on which representative of the isotopy class 
we choose. To make matters worse, aL) R 'y and a'U R 7 are 
not in the same isotopy class. However, these two ap- 
parent problems cancel each other out, as a result of the 
second constraint in l|58|) . (It's useful to draw a couple 
of pictures to sec this.) 

It is clear from the preceding considerations why the 
second constraint in (|58|l is necessary. The first con- 
straint in (|58l) is necessitated by consistency with the 
second one, as may be seen from fig. 

In defining the action of the Wilson loop operators, we 
ignored the fact that they are defined for parameterized 




FIG. 4: The second constraint, which assigns a —1 to the 
operation of cutting and rejoining of two parallel strands re- 
quires that a contractible loop change the wavefunction by 
-1. 



curves 7(5). This definition is manifestly invariant under 
orientation-preserving reparameterizations. Orientation- 
reversing reparameterizations simply conjugate the Wil- 
son loop operator: the integral is done in the reverse 
direction. For the case m = 2, the Wilson loop operators 
have real eigenvalues 1,-1, so they are also invariant un- 
der orientation-reversal. Thus, we can be lazy and treat 
the 7s as unoriented curves. 

The inner product on Hilbert space is determined 
by the condition that W [7], £[7] be unitary. This 
can be accomplished by inheriting the obvious inner 
product from pre-Hilbert space or by simply taking 

Thus, we have learned how to create a representation 
of the basic operator algebra l|49l) . I|5UII which defines 
J7(l)2 x U (1) 2 Chern-Simons theory in terms of simple op- 
erations on multi-curves, so long as we impose a carefully 
chosen constraint structure on the Hilbert space. Such 
a Hilber t spa ce necessarily realizes the doubled theory 
17(1)2 x U(l) 2 because it automatically admits a set of 
commuting operators L[y] as defined in the third line of 
(|60|l . Such a set of operators is a feature of the doubled 
theory, according to (|49(l . but not of the undoubled the- 
ory U(l)2- The undoubled theory suffers from a chiral 
anomaly, as we discuss in the penultimate section. Pic- 
tures of curves on surfaces cannot give rise to such an 
anomaly, so any theory with such a pictorial representa- 
tion must have no anomaly. 

Note that the Chern-Simons constraint (|34|l is already 
implemented in the pre-Hilbert space by the condition 
that states depend only on the isotopy class of a and not 
on a itself. This structure will be common to all of the 
topological field theories which we will dis cuss. The ad- 
ditional structure (|58|) is special to U(l)2 x U(l) 2 Chern- 
Simons theory, and it enables Ti. to realize the canonical 
commutation relations of the theory. In more complex 
topological field theories, these relations must be gener- 
alized. It will turn out that this can only be done in 
a rather restricted set of ways. In section IV, we will 
discuss this in detail. 

First, let us carry a little further our analysis of the 
1/(1)2 x U(l) 2 Chern-Simons theory. Let us consider the 
theory on the plane, but with a quasiparticle located at 
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the origin. 

S=^- I d 2 xdte»" p aiJl d u a p -^- [ d 2 x dte pvp c p d u c p 
An J 47r J 

+ J dt { Pl a (0,t) + p 2 c (0,t)) (61) 

We have coupled the gauge fields to a fixed, non- 
dynamical quasiparticle at the origin. The quasiparticle 
has charges p\2 under the two U(l) gauge symmetries. 
If pi = p% = 0, then there's no quasiparticle at the ori- 
gin. If one of them is non-zero then there is a non-trivial 
quasiparticle at the origin. These four possibilities are 
the four (topologically) distinct types of quasiparticles 
allowed in the theory, including the trivial particle. The 
Chern-Simons constraints now read: 

V x a = npi 

V X C = TT P2 (62) 

or, equivalently, 

W~[yo] = e™" 2 (63) 

for any curve 70 which encircles the origin. 

This can be incorporated into our representation of 
Hilbert space in terms of multi-curves if we enlarge our 
pre-Hilbert space to include multi-curves which termi- 
nate at the origin. In order to make this well-defined, 
we will also have to widen the origin into a finite-sized 
puncture in the plane. Let us also take the system to 
be on a finite disk, rather than the infinite plane. Thus, 
we must allow curves which are not closed but have two 
endpoints. The location of each such endpoint should be 
described by a quantum mechanical wavefunction ip(6) 
where 9 is an angular coordinate on the inner or outer 
circle. Different ips correspond to edge excitations, which 
we will discuss later. The endpoints of a curve may both 
be at the inner or outer boundary of the annulus or one 
may be at the inner boundary while the other is at the 
outer boundary. The former case corresponds to 'oscilla- 
tor modes' of the edge (which are, in general, gapped); 
the latter, to different sectors (or different Verma mod- 
ules) of the outer edge theory and different quasiparticle 
species at the origin. In the 1/(1)2 x U(l) 2 Chern-Simons 
theory which we are now discussing, the surgery relation 
tells us that we can reduce any configuration to one in 
which there is no more than a single curve connecting the 
inner circle to the outer one. For simplicity of depiction, 
we will choose preferred points at the inner and outer 
boundaries and insist that curves which connect the two 
boundaries terminate at these points, as shown in figure 
E] It doesn't matter which points we choose; different 
choices of preferred point are related to issues regarding 
boundary conditions and edge excitations, which we will 
discuss later. This is a particular choice of boundary con- 
dition. It is not natural physically, but it is convenient 



for now. Curves with two endpoints at the same bound- 
ary, or 'bigons', will be neglected for now and discussed 
in the context of edge excitations. 

From this construction, we see that the only allowed 
charges are p\2 = 0, 1. A Hilbert space representation in 
terms of multi-curves on the annulus does not exist for 
any other value of pi,2- This can be understood from a 
functional integral perspective by first noting from l|63|) 
that integer charges are only distinguished modulo 2. If 
pi is not an integer, then we can rescale to set p\ = 
1, thereby changing the coefficient of the Chern-Simons 
term. We really have a theory with a different (non- 
integer) coupling constant for one of the Chern-Simons 
gauge fields, and we shouldn't expect it to have a simple 
representation in terms of multi-curves. 




FIG. 5: Curves which terminate at the quasiparticle (i.e. the 
inner boundary) or the outer boundary of the annulus must do 
so at preferred points. The darker lines represent the bound- 
ary and the quasiparticle at the origin, while the lighter line 
represents a curve. 

Thus, we are essentially considering our system on an 
annulus. There are four basic pictures depicted in fig- 
ure EJ The first pair have been combined into two linear 
combinations with relative coefficient =Fl; the second pair 
have been combined into two possible linear combinations 
with relative coefficient ±i. All other one-dimensional 
submanifolds of the annulus (subject to the prescribed 
boundary conditions) can be obtained from these by ap- 
plying the operations b^qUO and )(— ► X. Thus, we 
can, as we did in the case of the torus, take a basis of 
states which vanish on the isotopy classes of all but one 
of the pictures in figure El It necessarily also vanishes 
on all isotopy classes which can be obtained from these 
three by the operations a — > a U Q and )(— > A given 
basis vector takes the value 1 on the isotopy class of one 
of the pictures and, therefore, takes the values ±1 on the 
isotopy classes of those pictures and those which can be 
obtained from it by the repeated use of the operations 
a — ► a U O an d )(—* We will call the four states 
corresponding to the four pictures in fig. El respectively, 
|0,0), |1,1), 1 1 , 0) , |0,1). Using @0), we can show that 
W + [7] has eigenvalues 1, —1, —1, 1 on these four states, 
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while W^~[7] has eigenvalues 1,-1,1,-1 if 7 encircles the 
origin. VF + [7] and IF - [7] have eigenvalue 1 on all four 
states if 7 does not encircle the origin. Hence, comparing 
with (|63[) . we see that our labels can also be interpreted 
as |pi, p 2 >- 




FIG. 6: The four species of quasiparticles (including the vac- 
uum as the trivial quasiparticle) in £7(1)2 x 17(1) 2 Chern- 
Simons theory. 

These particular linear combinations have been chosen 
so that they are spin eigenstates (by 'spin', we refer here 
to the eigenvalue under a rotation of the annulus, not to 
an internal quantum number, or, in more mathematical 
language, eigenstates under a Dehn twist). They also di- 
agonalize all two-particle braids. According to the spin- 
statistics connection, this is automatic for braids of iden- 
tical particles since we have taken spin eigenstates. We 
have furthermore chosen linear combinations such that a 
braid of two particles of two different particle species re- 
sults in merely the acquisition of a phase, i.e. we have di- 
agonalized two-particle braids. (The non-Abelian nature 
of the braid group can only be manifest when three or 
more particles are present.) As a result of these choices, 
our quasiparticles have the mathematical property of 
idempotence under 'stacking'. If we join annuli concen- 
trically (or stack the topologically equivalent cylinders) 
then (p 1 ,p 2 ) o (pi.p'a) = (pi, P2) 5 p lP { Sp 2 p> 2 , where o de- 
notes the stacking operation. 

Imagine taking two of these quasiparticles, with 
charges (pi,pi) and (p'x,^), very close together so that, 
viewed from a distance, they look like a single quasipar- 
ticle of charge (pi,P2,)- Following the graphical manip- 
ulations in figure [7| we find that p[ = (pi + p'j) mod2. 
This process is called fusion. 

We can also take one particle around another, as in 
figure El The effect of a counter-clockwise braid of two 
particles is the multiplication of the state by a phase 
e 2ie P1( ,[,p 2 p 2 _ e n(pip{-p 2 p' 2 ) _ Two identical particles can 
be exchanged counter-clockwise; the resulting phase is 




FIG. 7: The fusion of two particles. After the first step, the 
contractible loop can be shrunk directly or after applying the 
surgery operation. The same result is obtained either way 
(either 1 or 3 minus signs accrue). Strictly speaking, the 
particles are the superpositions of fig. [SJ so this picture must 
be superposed with 3 others; the result is the same. 

e l9 PiPi,P2P2 = e —(Pi~P-2) — e 1 f(pi-P2)^ An equivalent way 
of arriving at this result is via the spin-statistics theo- 
rem: we can imagine fusing the two particles first and 
then simply rotating the resulting particle by it. Under 
a rotation, a quasiparticle of charge (pj , P2) acquires a 
phase e mS — e~( Pl ~ P2 ) which is the same phase which 
we would obtain if we exchanged them first and then 
fused them. 

To summarize, we can se t up the Hilbert space not 
only of pure £7(1) 2 x U(l) 2 Chern-Simons theory, but 
also of the theory with static quasiparticles in a pictorial 
representation. 

The simplest generalizations of the preceeding are the 
other Abelian doubled Chern-Simons theories, U(l) m x 
U(l) m . These also have a pictorial representation, but 
it involves directed multi-curves and m-valent vertices at 
which they can terminate. We will not discuss them fur- 
ther here, but they are a straightforward generalization 
of the case m = 2. 



C. Z-2 Gauge Theory 

As we mentioned in the introduction, discrete gauge 
theories also have topological phases. Consider the sim- 
plest, gauge theory. There are two different ways of 
realizing such a theory. We could begin with a U(l) 
gauge theory with Maxwell action which is coupled to a 
charge-2 matter field. When this matter field condenses, 
the U(l) symmetry is broken to Zi- This construction 
can be done directly in the continuum. Alternatively, one 
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FIG. 8: The result of braiding two (1, 1) quasiparticles. (1, 1) 
quasiparticles are the type which are displayed in the lower 
part of fig. |S| with the plus sign. They involve a superposition 
of 2 pictures, hence two quasiparticles involve a superposition 
of four pictures. A clockwise interchange (second panel) leads, 
after surgery (third panel), to a phase tt/2. 

can work with Z2 gauge fields from the beginning. How- 
ever, one must, in such a case, work on a lattice. Let us 
follow the latter avenue. We consider a 2 + 1 dimensional 
space-time lattice on which there is an Ising gauge field 
degree of freedom a z = ±1 on each link of the lattice. 
We will label them by a lattice site, x, and a direction 
i = x,y,r so that there are three links associated with 
each site. The action is the sum over all plaquettes of 
the product of a z s around a plaquette: 

S = -K 2J Vz<Jz<J z cr z (64) 
plaq. 

To quantize this theory, it is useful to choose Coulomb 
gauge, cr z (x, t) — 1 for all x. In this gauge, the Hamilto- 
nian takes the form: 

H = -y"V x (x, i) - K a z a z a z a z (65) 

x - { spatial plaq. 

In Coulomb gauge, there are residual global symmetries 
generated by the operators 

G(x) = a x (x,x)a x (x,y)a x (x--k,x)a x (x.-y,y). (66) 

The extreme low-energy limit, in which this theory be- 
comes topological, is the K — > 00 limit. In this limit, 
GzVzVzOz — 1 for every spatial plaquette, recovering Ki- 
taev's 'toric code—. 



It is useful to define operators W[j] associated with 
closed curves 7 on the lattice: 

m= TL q (67) 

We also need operators Y[a] associated with closed 
curves on the dual lattice, i.e. closed curves which pass 
through the centers of a sequence of adjacent plaquettes. 

Y[a] = [] ( 68 ) 

x,ila 

The product is over all links which a intersects. L[ 7 ] is 
analogous to a Wilson loop operator while Y[j] creates 
a Dirac string. 

Let us consider the space of states which are annihi- 
lated by the Hamiltonian; this is the Hilbert space of 
the K — > 00 limit. When restricted to states within this 
Hilbert space, L[y] and Y[a] satisfy the operator algebra 

L[ 1 ]Y[a] = {-l) I ^Y[a]LU 

[L[ 1 ],L[a]] = [Y\ 1 ],Y[a]]=Q (69) 

Now, it is clear that such an operator algebra can be 
represented on a vector space which is very similar to 
the Hilbert space of [7(1)2 x ^(1)2 Chern-Simons theory: 

£M *[{«}] = *[{«}] 

Y[l] *[{«}] - *[{«U 7 }] (70) 

The notable difference is that the allowed states must 
now satisfy the constraints 

*[{«}] = *[{«uO}] 

*[{«}] = (71) 

Again, a is obtained from a by performing the surgery 
operation )(— > ^ on any part of a. 

If a is contractible, then Y[a] commutes with all other 
operators in the theory, so its effect on any wavefunc- 
tion should be multiplication by a scalar. If we take this 
scalar to be 1, then we have the first constraint above. 
The second constraint is necessary in order to realize the 
operator algebra (|70l) and is also required by consistency 
with the first. As a result of the second line of l|7ip. 
CVU7 in H70|) can be either cvUa.7 or uUlJ since they are 
equivalent in the low-energy Hilbert space. 

Again, we can characterize Hilbert space as the space 
of states annihilated by two projection operators, 

Ki = (lW)-|{«UO») (({«}|-{{aUO}|) 
P2,i = (]{<*}) - |{5}>) ((Ml - (72) 

As a shorthand, we will summarize such relations in the 
manner shown in figure 
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= l 



FIG. 9: As a shorthand, we will denote the constraints satis- 
fied by the physical Hilbert spac e as sh own. In the top panel, 
we have the relations of 17(1)2 x U(l) 2 Chern-Simons theory, 
while in the bottom panel, we have the relations of Z2 gauge 
theory. 



Again, the Hilbert space on the torus is four- 
dimensional and there are four quasiparticle species cor- 
responding to them. The corresponding pictures are the 
same as in Z7(l) 2 x U(l) 2 Chern-Simons theory, namely 
the four pictures in fig. However, as a result of 
the + sign in the second line of (|T1|) . the coefficients 
are different. =p in the top panel of fig. El becomes ±, 
while ±i in the second panel becomes simply ±. The 
fusion rule for quasiparticles (p\,p2) and (pi,/?^) is still 
pf = (pi + p'j) mod2. However, these quasiparticles have 
spins e 2 = e i7rpiP2 . A counter-clockwise braid of two 

quasiparticles yields a phase e l7T ( PlP ' 2+p ' lP2 ) . 

There is a natural generalization of this construction to 
discrete G gauge theories. In these theories, the variables 
are elements gi of the discrete group G on each link i of 
the lattice. The action is then of the form 



S = -K ^2 Tr (51525354) 
plaq. 



(73) 



We will not discuss these theories further here, but they, 
too, have a representation similar to the ones we have 
constructed, but with some extra features (e.g. directed 
curves) of the type present in the U(l) m x U(l) m Chern- 
Simons theories. 



D. Some Comments on Doubled Theories 

The two Hilbert spaces which we have just constructed 
are representation spaces for the underlying Wilson loop 
algebras of the U(V)2 x f7(l)a Chern-Simons theory, 

L[l]W + [i] = W/ + [ 7 ']L[ 7 ] 
L[-f)L[a) = L[a\L[j] 
W + [~f}W+[a} = {-iy^^W + [a}W + [~f} (74) 

(since W~ can be recovered from L and W + , they form 
a complete set) and the Z2 gauge theory, 

L[ 7 ]Y[a] = (-l)^y[a]i[ 7 ] 
L[ 7 }L[a] = L[a]L[j] 

Y[ 7 ]Y[a] = Y[a]Y[j\ (75) 



As we noted earlier, it is tempting to think of the L[y]s 
as 'number' operators and the VT + [7]s or Y[y]s as being 
somewhat similar to creation/annihilation operators. In- 
deed, L[y] essentially counts how much flux is enclosed 
by the curve 7, while VF + [7] and Y[y] can increase or 
decrease this flux. The only difference between the two 
theories is in the third lines of (|74|) and J7SJ: the V[7]s 
commute with each other while the IT + [7]s do not. 

This leads to rather significant differences between the 
two theories althoug h both have degeneracy 4 on the 
torus. The (7(1)2 X U(l)2 theory has in its spectrum 
of quasiparticles two species which are semions and anti- 
scmions. Z2 gauge theory has no quasiparticles with non- 
trivial self-statistics. The only non-trivial statistics are 
off-diagonal statistics between different particle types. 
Mathematically, these differences follow from the sim- 
ple fact that the composition algebra for Wilson loops 
on T 2 in the Z2 gauge theory is commutative and, there- 
fore, is equal to ©| =1 C; in the f/(l)2 x U(l)2 theory, it 
is is isomorphic to M2, the algebra of 2 x 2 matrices. 
(Any finite-dimensional C* algebra must be a direct sum 
matrix algebras* 8 , so these are the only possibilities.) 

In fac t, both of these theories are 'doubled' theories. 
[/(1)2 x /7(1)2 Chern-Simons theory, in the naive way: 
it is the tensor square of t/(l) 2 Chern-Simons theory. 
Z2 gauge theory is an example of a non-trivial double. 
From a formal perspective, it is the double of a theory 
with a single gauge field with vanishing Chern-Simons 
coefficient. Such a theory is gapless, and might represent 
a superfluid. 

We wil l procee d momentarily to a discussion of the 
SU(2)k x SU(2)k Chern-Simons theories, but it is worth 
pausing for a moment to note some of the general fea- 
tures of the theories which we have just constructed since 
these features will reappear in a slightly more compli- 
cated guise. We constructed our theories by finding a rep- 
resentation for the algebra of Wilson loop operators. We 
found that this could be done on a subspace of the vec- 
tor space of functionals of 1-manifolds, a subspace which 
was selected by requiring that allowed wavefunctionals 
assign the value d = ±1 to a contractible loop and be 
invariant under a surgery procedure on 1-manifolds. The 
doubled SU (2)& Chern-Simons theories theories have dif- 
ferent values of d and surgery procedures. 



SOLUTION OF DOUBLED SU(2) k 
CHERN-SIMONS THEORIES 



A. Wilson Loop Algebra 

We now turn to our main concern in this paper, the 
SU(2)k x SU(2)k Chern-Simons theories. Our aim is to 
construct the Hilbert spaces of these theories in a rep- 
resentation similar to the ones which we just used for 
L7(l) 2 x U(l)2 Chern-Simons theory and Z 2 gauge the- 
ory. 
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The action of SU(2)k Chern-Simons theory is 



In the spin- 1/2 representation, 



k f ( ,2 
— / tr [a A da H — a A a A a 
4tt ./ V 3 



a l T AB' a 2 



6 



(76) 



_- 27r cah a 6 

.2tt 

= Z T' 3 



[^d^bc - <5ac<5bd] (81) 



a.^ is a gauge field taking values in the Lie algebra su(2), 
a°T a , with index a — 1,2,3 running over the generators 
T 1 , T 2 , T 3 of su(2). In the second line, we have rewritten 
the action in the more compact language of differential 
forms. The integer, k, is the coupling constant. 

In constructing physical observables, we must exer- 
cise a little more care than in the Abelian case because 
the gauge fields at different points will not commute, 
[af l (x 1 )T a , al(x 2 )T b ] ^ 0. Thus, the exponential inte- 
gral must be path-ordered 



U[j] = Ve^ TC - dl 

oo 

E 



2rr />s± 

/ I ds\ I ds2 ■ ■ 



ds n 



7(*i)-a& ( 7 (*i)) T a > . . . j(s n )- a ^ (j(s n )) T a " 



(77) 



where 7(s), s S [0, 27r] is an arbitrary parameterization 
of the curve 7. This quantity is an SU(2) matrix, which 
transforms in the adjoint representation of SU (2) at the 
starting point 7(0). To get a gauge-invariant quantity, 
we must take the trace: 



We must now be a little more careful about issues of 
gauge-invariance. The above commutation relations hold 
in a field theory which contains 'too many' degrees of 
freedom, most of which are pure gauge. Continuing in 
parallel with our discussion of the Abelian case, we will 
work entirely with gauge-invariant Wilson loop opera- 
tors. In so doing, we are eliminating the pure gauge 
degrees of freedom. In the Abelian case, this is trivial 
since the gauge transformation is simply a — » a — df. In 
the non- Abelian case, it is a — > gag -1 — dg so there 
is a non-trivial Jacobian which results upon eliminating 
the gauge degrees of freedom. As a result of this Jaco- 
bian, the commutation relations of the reduced the ory 
are modified by the shift k — ► k + 2, as shown in ref. |2Jj. 
(The alternative possibility is to work with the full set of 
degrees of freedom of the theory and then require that 
the Hilbert space contain only gauge-invariant states. In 
such a case, the inner product will be non-trivial and will 
lead to this 'quantum correction'.) We will see how this 
shift arises from a different perspective in section IVT1 so 
we defer a discussion until then. 

We would like to work in this reduced theory, contain- 
ing only gauge-invariant degrees of freedom. Hence, we 
make this shift when we compute the commutator of two 
Wilson loop operators, which is: 



W[ 1 ]=tr{U[ 1 ]) 



(78) 



[W[7],W[Y]] = 



2 sin 



2(fc + 2) 



Y, (WfroiY] - W[ 7 'o i7 ]) (82) 



As in the £7(1)2 x f (1)2 case, the Wilson loop op- 
erators in this theory are, strictly speaking, defined for 
parameterized curves 7(5). It is manifestly invariant un- 
der orientation-preserving reparameterizations. However 
unlike in the U(l)2 x ^(1)2 case, it is not quite true here 
that W[y] is real, so the orientation of 7 cannot be treated 
so cavalierly. Fortunately, the spin- 1/2 representation of 
SU(2) is pseudo-real, meaning that it is equal to its con- 
jugate representation upon multiplication on the left and 
right by the antisymmetric tensor e a p. These factors dis- 
appear upon taking the trace, so we can, again, ignore 
the orientation of 7, i.e. if we orient 7 and denote by 7 -1 
the same curve followed in the opposite direction then 



Wl-r 1 } = tr (£/[ 7 ]t) = tr (eC/[ 7 ]e) 



tr(C/[ 7 ]) 
W[y] (79) 



Proceeding in parallel with the Abelian case, we derive 
the equal-time commutations relations from the temporal 
gauge form of the action (|76|l 



The summation is over all intersections i of the curves 7 
and 7'. In this expression, 70^7' is the curve obtained by 
starting at intersection i, following 7' until it returns to 
i, and then taking 7. In so doing, we temporarily intro- 
duce orientations for 7 and 7' so that their tangent vec- 
tors form a right-handed dyad at their intersection point. 
Alternatively, we could simply say that as the intersec- 
tion is approached along 7, we turn to the left to join to 
7'. This construction is closely related to the Goldman 
bracked. Strictly speaking, we should have a factor of 
(7,0:), multiplying 70^7', where (7, a)i is the sign of the 
i th intersection between 7 and a. However, since the ori- 
entation is unimportant for closed SU{2) Wilson loops, 
as we just noted, we can drop this factor. However, for 
loops terminating at boundaries or for other Lie groups, 
one must retain this factor. 

In the doubled theory, we have two mutually commut- 
ing sets of Wilson loop operators, W±[7], which satisfy 
the algebra 

(83) 



«T: a 2 



k 



(80) 



[W + [j],W + [j']} = 

(A-A-^YiW+boiY. 



18 



[W_[ 7 ],W-[Y]] = (84) 
- {A - A- 1 ) J2 (W- boa'} ~ W- [ 7 'o l7 ]) 

[W+H.WlfY]] = 1 (85) 

where we have introduced A — i exp (iri/2(k + 2)). Note 
that if 7 and 7' have multiple intersections then 70,7' 
will be a self-intersecting curve. 

To construct a Hilbert space which furnishes a rep- 
resentation of this algebra, we begin with a structure 
similar to that of the Hilbert spaces which we derived 
earlier: it is a subspace of the vector space of functionals 
of isotopy classes of loops which satisfy 

*[{aUO}] - <**[{"}] (86) 

where d is a constant to be determined. The second 
condition in l|58l) is replaced by a more complicated one 
which we will construct later. For the following discus- 
sion, it is useful to introduce the notion of e?-isotopy. Two 
multi-curves are related by d-isotopy if one can be de- 
formed into the other by a combination of isotopy and 
the elimination of contractible loops. The latter oper- 
ation results in a factor of d multiplying the wavefunc- 
tional. 

We will find it notationally convenient to define '5 [{a}] 
for intersecting loops so long as intersections are resolved 
as over- or under-crossings. In other words, we allow the 
1-manifold a to no longer be embeddable in a surface. We 
will relate in the following way evaluated on a with 
crossings to its values on smooth 1-manifolds embedded 
in a surface. We define: 

= A9[{a/}] + A-^Ha"}] (87) 

where a' and a" are the two ways of resolving the in- 
tersection, as depicted in figure ITUl This relation can be 
applied repeatedly to remove all crossings in a, We could 
have implemented this decomposition into resolutions in 
the definitions of W± [7] , but it is a little simpler this way. 



that 7 always crosses over a in 0*7; in 7*a, a always 
crosses over 7. 

We also need to define W+[{7}] for self-intersecting 
curves such as 7047'. A natural definition is to start at 
a given point on the curve, follow the curve in a given 
direction, and rule that 'later' sections of the curve al- 
ways cross over 'earlier' sections. In other words, we 
parameterize the curve as j(t), with t S [0, 2tt] and 
7 (0) = 7 (2tt). We say that if j(t) = 7(f), then >y(t') 
crosses over j(t) if t' > t. However, this depends on both 
the starting point and the direction. Hence, we aver- 
age over all possible starting points and both directions. 
Thus, W+ [{7}] for self-intersecting 7 is a normalized sum 
nEm^+K)™}] where the j m have intersections speci- 
fied as over-crossings or under-crossings. Note that this 
is not a sum over all possible ways of choosing the in- 
tersections to be over- or under-crossings. If there are 
n intersections, then there are n different 7™ which can 
result in this way, not 2". 

Before verifying that the operators defined above sat- 
isfy the desired commutation relations, let us first make 
sure that they are well-defined. Since \f[{a}] only de- 
pends on the isotopy class of a, W±[{j}] should only 
depend on the isotopy class of 7. Thus, 0*7 must be 
invariant under a continuous deformation of 7 into 7' 
which has two new intersections with a, as shown in fig- 
ure ^2 Applying the definitions (|88|l . we see that 0*7 is 
invariant under isotopy moves of a and 7 if: 




FIG. 11: a*7 must be denned so that 0*7 = 0*7' 





a 



a 



a 



FIG. 10: a' and a" are the two ways of resolving the crossing 
in a 

The action of the Wilson loop operators is given by 

^+[{7}] *[W] = *[{a*7}] 
VF_[{ 7 }j*[{a}j = *[{7*a}] (88) 

The operation * is defined as follows: if a and 7 do not 
intersect then it is simply the union a U 7. However, all 
intersections between a and 7 are resolved by specifying 



The isotopy invariance of W± [{7}]^ [{a}] follows from 
its close relation to the Kauffman bracked of the 1- 
manifold obtained by overlaying 7 on a (or the reverse, 
in the case of W- ) . The Kauffman bracket is defined for 
multi-curves in R 3 by projecting them to the plane, but 
keeping track of over-crossing and undercrossings. The 
crossings are resolved according to the rule l|87|l and all 
unknotted loops are accorded a factor of d. The only dif- 
ferent feature in our Hilbert space is that the loops are 
actually embedded in some surface so that some unknot- 
ted loops are not contractible in that surface. 

We now verify that W±[{7}] as defined in obey 
the commutation relations (|83|) . First, note that ^+[{7}] 
and W- [{7'}] commute trivially because 

IF+ [{7}] WM{ 7 '}] = *[{ 7 '*«*7}] 

= IF_[{ 7 '}]VF + [{ 7 }]*[{a}] (90) 
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Now consider the commutation relation between 
^+[{7}] and ^+[{7'}] (the situation for ^[{7}] and 
^-[{7'}] i s so similar that we need only discuss W + ). If 
there are no intersections between 7 and 7', then both 
the left- and right-hand-sides of l|83|l vanish. Suppose 7 
and 7' have a single intersection k. The left-hand-side of 
(|83[1 acting on a state ^[{a}] is 

LHS = *[{a*7'*7}] - *[{a*7* 7 '}] 

= A*[{a*( 7 'o fc7 )}] + A- 1 ^[{a*( 7 o fc7 ')}] 

-AV[{a * (7°fc7')}] - A~ l ^[{a * (7'°fc7)}] 
= (A-A- 1 )(*[{a*( 7 / o fc7 )}]- 

*[{a*(7o fc7 ')}]) 
= RHS (91) 

Consider now the case in which 7 and 7' have two in- 
tersections. In order to analyze this, it is useful to sym- 
bolically denote the two resolutions of the first crossing 
by x\, yi and the two resolutions of the second crossing 
by x 2 , 2/2- Then the left-hand-side of (|83|) is: 

LHS = {Ax l +A- l y l ){Ax2+A- 1 y 2 ) 
- (A yi + A- x Xl ) (Ay 2 + A~ l x 2 ) 
= (A 2 - A- 2 ) ( Xl x 2 - ym ) (92) 

Meanwhile, the right-hand-side is: 

RHS = (A- A- 1 ) (x 1 (A + A- 1 )(x 2 +y 2 )/2 

- Vl (A + A- 1 ) (x 2 +y 2 )/2^j 
= (A 2 - A- 2 ) { Xl x 2 - ym ) (93) 

so (I83|l is satisfied. 

If 7 crosses over 7' at both intersections (or vice versa), 
then we can deform either one so that there is no intersec- 
tion. Thus the equality which we have shown is trivial: 
both sides of the equation vanish because the correspond- 
ing Wilson loop operators commute, as we discussed ear- 
lier. If 7 crosses over 7' at one intersection and under 
it at the other, then the commutator will be non-trivial, 
but still satisfies 183(1. as we have just seen. The general 
case of arbitrary 7, 7' is similar. 

As in the Abelian case, we have found an infinite- 
dimensional vector space on which we can represent the 
commutator algebra of our theory. As in that case, 
we must truncate this vector space because the classi- 
cal phase space has finite volume. Consider the torus. 
The holonomy E/+[ 7 i] of the a+ gauge field about the 
meridian, 71, of the torus will be some SU(2) rotation. 
Since the homotopy group of the torus is Abelian, the 
holonomy about the longitude of the torus i7+ [72] must 
commute with [/ + [ 7 i], i.e. [t/+[ 7 i], [/+[ 7 2]] = 0. Hence, 
we must specify two SU(2) rotations about the same axis. 
The direction of this axis is not invariant under an SU (2) 
gauge transformation, so we need only specify two angles. 
The same is true for the a_ gauge field. Hence, the phase 
space of the theory is the product of two tori, T 2 x T 2 . 



Since the phase space of the classical theory has finite 
volume, its Hilbert space is finite-dimensional. In order 
to truncate our Hilbert space to a finite-dimensional one, 
we must specify surgery relations which must be obeyed 
by states in Hilbert space, analogous to the surgery re- 
lation of the d = — 1 theory shown in the lower panel of 
fig. For the level k theory, there is only one possible 
relation which is consistent with the corresponding value 
of d, the amplitude associated with a contractible loop. 
The value of d so strictly constrains the Hilbert space of 
our theory that we have no freedom at all in our choice of 
a surgery relation - the analogue of )( = — If we in- 
troduce no constraint, the Hilbert space on the torus (or 
any higher genus surface) will be infinite-dimensional. If 
we introduce a constraint which is too severe for the given 
d, such as ) ( = for d = \/2, then the Hilbert space (on 
any surface) will be zero-dimensional since there won't be 
any wavefunctions which satisfy both constraints. For a 
given d, there is a unique constraint which is just right 
- neither too trivial, nor too severe - so that there are 
finite-dimensional Hilbert spaces associated with differ- 
ent surfaces. In fact, for arbitrary d, there are no such 
non-trivial constraints, so the Hilbert space must be ei- 
ther infinite-dimensional or zero-dimensional. Only for - 
surprise, surprise - the very same sequence of ds which 
we have found in this section in 1891 d = 2 cos(7r/(fc + 2)), 
are there non-trivial constraints which lead to finite- 
dimensional Hilbert spaces. In section IV El we will find 
these constraints. 

As we will see in that section, such surgery relations 
depend strongly on d. We found above that d is de- 
termined by the condition that states and operators be 
invariant under isotopy which, in turn, follows from the 
Chern-Simons constraint which requires that the connec- 
tion be flat. However, it is useful to consider another 
perspective on how d is determined in SU(2)f- x SU(2)k 
Chern-Simons theory. In order to do this, it will be use- 
ful to have an inner product on our Hilbert space. We 
will initially define this inner product on our infinite- 
dimensional Hilbert spaces, but we will later show that its 
restriction to the finite-dimensional Hilbert spaces con- 
structed in section lYEl is the correct inner product of our 
theory. 

B. Inner Product 

The infinite-dimensional pre-Hilbert spaces of the pre- 
vious section have a natural inner product: 

<^}|M> - { I otWise ( 94 ) 

where = denotes equivalence under isotopy. The pleas- 
ant surprise about this inner product is that W± [7] are 
Hermitian with respect to it: 

{{a'}\W+[y}\{a}) = ({«'}({« * 7 }) 

= «M|{a'* 7 })r 
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= (({a}\W + W\W})Y (95) 

The second equality can be understood by noting that the 
inner product vanishes unless a' = aU#7 or a' = aUi7- 
However, a' = aU^7 implies that a = q;'Ul7, and vice 
versa. The equality follows since one is accompanied by 
a factor of A; the other, a factor of A^ 1 . 

Classically, W±[7] are real because the spin-1/2 rep- 
resentation of £{7(2) is pseudoreal, as we noted in the 
previous section. Thus, it is natural to demand that 
they be Hermitian operators in the quantum theory and, 
indeed, we could take this requirement as the defining 
condition of our inner product. Fortunately, as we have 
just seen, this leads to precisely the same inner product 
as the natural one on pre-Hilbert space. Eventually, we 
will restrict this inner product to the finite-dimensional 
subspace which will form the Hilbert space of our theory. 

C. Accidental Symmetry and Some Coincidences 

Before giving any detailed calculations, some low-level 
'coincidences' should be brought to light so that they do 
not cause confusion later. These are essentially of the 
same form as the 'accidental' SU (2) symmetry of a free 
boson at the self-dual compactification radius R = l/v2 
which underlies the Abclian bosonization of an SU(2) 
doublet of fermions. A free chiral boson <p has the con- 
served current idip. The existence of this current en- 
dows it with a U(l) Kac-Moody algebra. The theory has 
central charge c = 1 (with respect to the Virasoro alge- 
bra) for any compactification radius R, tp = tp + 2ttR, so 
the description in terms of a U(l) Kac-Moody algebra 
is perfectly acceptable. However, at radius R = 1/V2, 

there are additional dimension-1 fields e ±lv ^ - cur- 
rents - which are allowed by the angular identification. 
These three currents form an SU{2)\ Kac-Moody alge- 
bra. Thus, at this special radius, we can describe the 
theory equally well in terms of an SU{2)\ Kac-Moody 
algebra, which also has central charge 1 with respect to 
its enveloping Virasoro algebra. 

Let us now turn to the Abelian theories discussed in 
previous sections, (A) U(l)2 x 1/(1)2 and (B) Z 2 gauge 
theory. They are 'doubled' SU (2)i theories in a sense 
which we now describe. Consider the action of a Wilson 
loop operator in doubled SU{2) Chern-Simons theory. 

W+ [{7}] *[{«}] = *[{a*7}] (96) 



Consider an intersection between a and 7. Using the 
prescription l|87|l for resolving overcrossings, 

W+IM] *[{a}] = A%{a\J Ll }] + A-^aUuf}] (97) 

If we now take A = e 7 ™/ 6 , which is almost but not quite 
what we expect for k = 1, and use the surgery relation 
)(= — (we will not show that this is the correct relation 
for the k = 1 theory until section lV El but let us go ahead 
and use this relation nevertheless) then we find that this 
is simply 

W+ [{7}] *[{a}] = 2isin(7r/6)*[{aU L7 }] 

= i^[{aU Ll }} (98) 

Furthermore, d = —A 2 — A ~ 2 = — 1. Thus, this theory 
is equivalent to the U(l)2 X U(l)2 Chern-Simons theory. 
If, on the other hand, we take A = i e 7 ™/ 6 , as expected 
for k = 1, then we find d = 1, which leads to Z2 gauge 
theory. Thus, doubled SU (2)i is Z 2 gauge theory, while a 
slightly modified version is U(l)2 x £7(1)2. Chern-Simons 
theory. 

In section IV El we will show that these theories are 
so tightly constrained by combinatorial relations, so that 
self-consistency essentially specifies the entire structure. 
In the combinatorial worloj 5 ^* 5 ^, at level k the SU(2)k the- 
ory is fully specified once a primitive 4?* th root of unity 
A, with r = k + 2, is given. (If A is only a primitive 2r th 
root, the modular S-matrix is singular, but the Drin- 
feld double - as opposed to the mere tensor square of 
the singular theory - does have a non-singular modular 
S'-matrix.) Our theory A is the doubled SU "(2)i theory 
for A = e 2 W 12 w hile theory B is the Drinfeld doubled 
SU(2) 1 theory for A = ie 2 ™/ 12 . Both of these have level 
1 as 577(2) theories and, according to the Kauffman re- 
lation d = -A 2 — A~ 2 have d = — 1 and d = +1. So, 
although we introduced these theories through their re- 
lation to the Abelian groups U(l) and Z2, they are also 
low-level SU(2) theories. 

For the combinatorially defined SU(2) theories (and 
their doubles) to be unitary, there is a strong restriction 
on A. In fact, for k > 1, A must be chosen as A = 
±ie 27Tl / 4r , r = k + 2, in order that the intrinsic inner 
product, even on the closed surface Y of genus 2, be 
positive 53 . The unitarity of our first example (theory A 
of section IV C|l is a bit of an exception as seen in the 
table of undoubled SU(2) theories below. 





k even 


k = 1 


k > 3 odd 


A = e 2 ™/4(fc+2) 


non-unitary 


unitary (d=-l) 


non-unitary 




S'-matrix non-sing. 


S-matrix non-sing. 


S-matrix non-sing. 




unitary 


unitary 


unitary 


^4 _ l e 27T4/4(fc+2) 


S-matrix non-sing. 


S-matrix singular 


S-matrix singular 






but non-sing. S matrix for 


but S non-sing, if restricted 






doubled theory=Z 2 gauge theory 


to integer spins 
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D. Contractible Wilson Loops: the Value of d 

In section IV Al we saw that isotopy invariance fixed 

d = 2 cos in the level k theory. This is somewhat 

surprising since d is the eigenvalue of W± [O] m any state 
in the theory. Since the gauge fields a± are flat, we would 
expect their holonomies to be trivial, t/ipy] = I2 and, 
therefore, W±[7] = 2. Why is d reduced from its naively 
expected value? 

In our discussion of the Abelian theories, we were able 
to get away with a rather fast and loose treatment of Wil- 
son loop operators. In general, these operators need to be 
regularized. To see why, consider the expectation value 
of a product of Wilson loops in the full 2 + 1-dimensional 
theory: 

(W[n] . . . W[ ln ]) = J VaW^x] . . . W[j„] e Scs (99) 

In an Abelian theory, this is equal to 
ex P(lr Ey L (lhlj)) wh ere £(7^7j) is the linking 
number of 7$ and 7,-. The problem is caused by the 
i = j terms. The self-linking number is not well-defined 
without some kind of regularization, e.g. point-splitting. 
In mathematical terms, this is called a framing of the 
curve 7i. One thickens the curve 74 into a ribbon and 
then computes the linking number of the curves at the 
two ends of the ribbon. Clearly, this is not unique, 
since the ribbon can twist an arbitrary number of times, 
but once a framing has been chosen, a well-defined 
calculation can be done and the result for different 
choices of framing can be related to each other. In more 
physical terms, the amplitude for such a process depends 
not only on how the different anyons wind around each 
other but also - since each has fractional spin - on how 
each particle rotates during the process. This extra 
information - which is equivalent to the framing - must 
be specified in order to have a well-defined process. 

In an Abelian theory, we can always choose a framing 
so that an unknotted contractible loop has self-linking 
number zero. Thus, it is possible to ignore this subtlety. 
In our construction of [7(1)2 X ^(1) 2 Chern-Simons the- 
ory above, we actually took an unknotted loop to have 
self-linking number 1 since W+[7] x I'[{a}] = ^[{aU^}] = 
— ^[{a}] = e 7rl l 5 , [{o;}]. In a non-Abelian theory, it is 
not possible to choose a framing so that the value of a 
Wilson loop is always unity if the loop is an unknotted 
contractible loop. It will, however, always be some con- 
stant, which we have called d: 

* [{«UO}] = rf* [{a}} (100) 

Let us compute d in SU(2)k Chern-Simons theory, fol- 
lowing the arguments of WittenSi. We will find that it 
is real, so its value will clearly be the same in the tensor 



square of the theory, which is our ultimate interest. Con- 
sider Chern-Simons theory on D 2 x R - space is a disk 
D 2 and R is the time direction - and a curve 7 E D 2 x R. 

It may seem that the functional integral H99|l. even 
when evaluated for a single Wilson loop W[y] is not quite 
the same thing as the operator W[y] which is the starting 
point for our canonical quantization procedure because a 
curve 7 in the functional integral need not lie on a con- 
stant time spatial slice. However, the distinction is illu- 
sory if the curve is unknotted because the Chern-Simons 
action is independent of the spacetime metric. Thus we 
can foliate spacetime into 'spatial slices' in any way that 
we like, in particular so that our unknotted curve 7 lies on 
such a slice. This can be said slightly differently by not- 
ing that the value of (|99|l for an unknotted contractible 
loop does not depend on whether or not it lies in a sin- 
gle spatial plane, it only depends on its topological class. 
Hence, a computation of W[j\ using the functional in- 
tegral should be the same as the result obtained from 
canonical quantization. 

For convenience, let us assume that 7 does not lie in 
a single spatial slice D 2 and that it intersects any spa- 
tial slice in either 2 or zero points (except for the two 
spatial slices which are tangent to 7). We define another 
topologically trivial curve 7' which is simply a copy of 
7 translated spatially. Since the functional integral only 
depends on the topological class of the curves 7, 7', we 
can take them far apart so that (W^W^']) decouples 
into (W[i\) {W[Y]). Thus, d 2 = (W[-y]W[y']) . 

Let us divide D 2 x R into two halves, D 2 x (— 00, 0] 
and D 2 x [0, 00) such that 7, 7' each intersect the spatial 
slice t = at two points. We will call these points xi, x 2 
and x' 1; x 2 , as depicted in figure El Each is divided by 
this slice into two arcs, which we call 7_, 7+ and 
Then we can define a state in the t = Hilbert space by 
performing the functional integral 

il>[A(x)]= I T>a(x,t) W[t_] W[i_] x 

Ja(x,0)=A(x) 

J^dtfd^Ccs (101) 

We can also define the state 
x[A(x)} = f Va{x,t)W[ 1+ ]W[i + ]x 

J a(x,0)=A(x.) 

e f Q °°dtfd 2 x Cos ( 1Q2 ) 

The inner product of these two states is the functional 
integral which we would like to compute: 

d 2 = (x|V»> 

= J Va{x,t) W[y] W[i] e t^ dt -f d2xCcs (103) 
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Consider, now, the state obtained by deforming 7_ , y_ 
in order to perform a counterclockwise exchange of X2 
and x'[. This state, which we will call B\tp) is depicted 
in!12l Consider, as well, £? _1 |?/>), obtained by deforming 
7_, in order to perform a clockwise exchange of X2 
and x' x . From the figure, we see that 

( X \B\ip) = d (104) 
(XIB- 1 ^) = d (105) 

In section IV Fl we will show that the four-quasiparticle 
Hilbert space is two-dimensional. Thus, B has two eigen- 
values, Ai, A2, so that 

B - {\i + \ 2 ) + \ 1 \ 2 B- 1 = (106) 

which, in turn, implies that 

d- (Ax + \ 2 )d 2 + X l \ 2 d = (107) 

so that 



In section I VII we will calculate these eigenvalues and 
show that 

d = 2cos {kT2) < 109 > 

E. Truncation of Hilbert Space: Jones- Wenzl 
Projectors 

In this section, we take up the issue of truncating the 
pre-Hilbert spaces of section IV AI to finite-dimensional 



ones. We will do this by finding the analogues of the 
second equations of l|58l) and (|71|l which reduced the re- 
spective Hilbert spaces of U(l) 2 x U(l) 2 Chern-Simons 
theory and Z 2 gauge theory. 

We begin by considering the structure of such con- 
straints in general. Suppose that we wish to impose a 
relation involving n strands of a given 1-manifold (i.e. n 
segments of the 1-manifold which are 'close together'). 
We assume that there are no relations involving fewer 
than n strands; if there were, we could always use it to 
reduce a set of n strands to fewer than n strands, and 
the n-strand relation would be superfluous at best and 
incompatible at worst. Now consider our n-strand rela- 
tion. If we were to connect the endpoints of two of the 
strands, then we would have ann-1 strand relation. By 
assumption, this is impossible. Hence, the putative n — 1 
relation must actually vanish identically. The same must 
be true for any other way of connecting two strands to 
yield an n — 1 strand relation. Needless to say, any n — 2, 
?i — 3, ... strand relations obtained in such a way must 
also vanish identically. This is a severe condition on our 
n strand relation. 

In order to construct relations which satisfy this condi- 
tion, it is useful to introduce the Temperley-Lieb algebra, 
TL n . This algebra is most simply described in pictorial 
terms. The Temperley-Lieb algebra on n curves is made 
up of all planar diagrams without intersections in which 
n curves enter at the bottom and exit at the top. Some 
elements of the Temperley-Lieb algebra are depicted in 
the top panel of fig. The algebra is generated by 

1, ei, e-2, . . . , e n _i, depicted in the bottom panel of fig. 
IT31 In the identity element of the Temperley-Lieb alge- 
bra, all curves go straight through. All other elements of 




FIG. 13: In the top panel, three elements of TL 5 are de- 
picted. In the bottom panel, the generators 1, ei, ea, - . - , e n _i 
are shown. 

the Temperley-Lieb algebra involve 'turnarounds' which 
join two incoming or two outgoing curves. The multipli- 
cation operation is simply stacking two pictures. Addi- 
tion is simply formal linear superposition, i.e. pictures 
are multiplied by arbitrary complex numbers and super- 
posed. The generators a satisfy the defining relations 

e\ = de, (110) 
ei ej — ej ei for \i — j\>2 (HI) 
e-i&i±\e-i = ei (112) 

Our desired n-curve relation is some element of the 
Temperley-Lieb algebra. Call it P n . The condition that 
all n — 1, n — 2, n — 3, ... curve relations which can be 
derived from it must vanish identically is the statement 
that P n annihilates every generator except the identity. 
The P n are known as the Jones- Wenzl projectors. 

Such projection operators are unique, up to an over- 
all scalar factor. To see this, imagine that there were 
two such operators P n and P' n . Since P n and P' n are 
themselves elements of the Temperley-Lieb algebra, they 
can be written as P n = 1 + / and P' n = 1 + g (we 
are free to choose a scalar factor, so we set the coeffi- 
cient of the identity to one in both expressions). Then 
P n P' n = Pn(l + 9) = Pn- However, it is also equal to 
P n P' n = (1 + f)P' n = P' n . Hence, P n = P' n . 

We can construct the P n s recursively. In order to do 
this, we first define the numbers A„, which are the traces 
of the P n s: we join every curve coming out of the top with 
its partner at the bottom, as shown in figure lT^I and eval- 
uate the resulting diagram by assigning a factor of d for 
each closed loop. Now, P 2 can be found by inspection, 
P2 =)( — 4 ^ (look familiar?). Let us suppose that we 
know P n -i, Pn—2i Pn-Si- ■ •■> P%- Then, connect one of the 
curves entering P n _i with the corresponding curve leav- 
ing the top. The result annihilates all turnarounds on 
n — 2 curves, so it is proportional to P n —i- By comparing 
their traces, we see that the constant of proportionality 
is simply the ratios of the traces A n _2/A„_i. Now, con- 
sider the element of TL n depicted in figure El Clearly, 
the n — 2 turnarounds which are entirely on the first n—1 
curves are annihilated. Consider a turnaround on the last 
two curves. As we see from the picture if we note that 
P n -\Pn-2 = Pn-i, it, too is annihilated. Thus, the oper- 
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FIG. 14: The trace of an arbitrary element (the shaded re- 
gion) of the Temperley-Lieb algebra. 



ator shown in figure ^| is, indeed, the desired projection 
operator. 




FIG. 15: In the top panel, the Jones-Wenzl projector P n 
is written in terms of P n -i- In the bottom two panels, it 
is shown how P„ annihilates a turnaround on the last two 
curves. 

Now let us apply this result to the problem of finding 
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an n-curve relation for a topological field theory with a 
pre-Hilbert space of states [{a}]. As we noted above, 
if we connect any adjacent endpoints out of the set of 2n 
cndpoints of n curves, then the ensuing relation must 
vanish. If we divide the 2n endpoints into n incom- 
ing ones and n outgoing ones then it must not only be 
true that any turnaround attached to the top or the 
bottom or bottom must be annihilated, but also that 
connecting the first top endpoint to first bottom one 
or the last top endpoint to the last bottom one yields 
zero. Imposing the latter requirement on P n implies 
that dP n -i ~ -Pn-i = or, simply A„_ 2 (d) = 

dA n -i(d). We have written the d dependence explic- 
itly in order to emphasize that this condition restricts 
d. Only those ds which satisfy this relation can lead 
to a consistent n-curve relation. By taking the trace of 
our equation for P n , we see that A„ = dA n _i — A„_2- 
Hence, our condition is simply that the trace of P n van- 
ishes. If we pick a d such that A n vanishes but A„_i, 
A„_2, etc. don't vanish (which will, in general, be pos- 
sible) then we have the desired n-curve relation. 

The equation A„ = d A n _i — A„_2 is the recursion re- 
lation for the Chebyshev polynomials. Only at the roots 
of the n th Chebyshev polynomial is there a consistent n- 
curve relation. Fortunately, the non-trivial roots of the 
(k + l) th Chebyshev polynomial are precisely the num- 
bers dk = ±2 cos (t^2^ which arise in SU(2) k x SU(2) k 
Chern-Simons gauge theory. 

Thus, there is only one possible relation which we can 
impose in the level- fc theory, and it involves k + 1 strands. 
We are 'lucky' that we could even impose one. For other 
values of d, there isn't even one such consistent relation. 
Let us write our k + 1-strand relation as 

l + ci/i + c 2 / 2 + ... + cJ g = Q (113) 

where 1, /i, /2, . . . , f q are elements of TLk+i- Consider 
a 1-manifold a in which k + 1 parallel strands have been 
brought close together. (In order to state this precisely, 
we may need to temporarily introduce a metric, define 
'close', and then show that, as a result of isotopy in- 
variancc, the final answer is independent of the metric.) 
Then let us define fj ■ a as the 1-manifold which is ob- 
tained by replacing the k+1 parallel strands by the ele- 
ment fj of the Temperley-Lieb algebra TL^+i. 

Using this notation, we can write the Hilbert space of 
SU(2)k x SU(2)k Chern-Simons gauge theory as the vec- 
tor space of functional of isotopy classes of loops which 
satisfy 

*[{aUO}] ~ dk *[{«}] - 
*[{«}] + = (H4) 

3 

where dk = 2cos(^-^^j, and v?[.] is defined for multi- 
curves with over- and under-crossings according to eq. 

E3 



Applying the second relation, we see that we can re- 
late \f r [.] evaluated on an arbitrary isotopy class to 
evaluated on isotopy classes with winding numbers less 
than or equal to k. Thus, the Hilbert space on the torus 
is (k + l) 2 -dimensional. 

Let us consider, for the sake of concreteness, the 
case k — 2. From the above discussion, we see that 
the Jones- Wenzl projector P3 can be constructed from 

P2 =)( - 2 X ' usin S A 2 = d2 ~ !> A i = d, and d = \/2. 
It is displayed in figure There are five terms in P 3 , 
and the corresponding rel ation sat isfied by states in the 
Hilbert space of SU(2) 2 x SU {2)2 Chern-Simons theory 
is shown in figure [TBI 

- /2¥[ 1 - /2¥[ 1 
+ V[ 1 + V\ 1 = 

FIG. 16: Th e condit ion imposed on states in the Hilbert space 
of SU "(2)2 x 5(7(2)2 Chern-Simons theory by the Jones- Wenzl 
projector P3. 

Thus, the Hilbert space on the torus is given by states 
^(nim)[{ a }] which vanish on all multi-curves a except 
those which are d-isotopic to the multi-curves (n, m) 
depicted in figure 1171 For these multi-curves, the d- 
relation and isotopy can be used to determine their 
value relative to the value of ^ ( n , m ) [•] evaluated on some 
fixed state in the d-isotopy class. One might have ex- 
pected the states to correspond to the set of (n, rn) with 
n, m = 0, 1, 2, but instead of (2, 2), we have (—1, 1) since 

*(2,2) = 2 V2*(0,0) _ *(2,0) _ *(0,2)- 




FIG. 17: The 9 ground states of the k — 2 theory on the torus. 
Opposite sides of the square are identified, as indicated by the 
arrows. 



F. Quasiparticles in SU(2)k X SU(2) k Chern-Simons 
Gauge theory 

Let us now consider quasiparticles and their braiding 
statistics. As in the Abelian cases which we consid- 
ered earlier, quasiparticles are modeled as interior bound- 
aries or punctures at which curves can terminate. Thus, 
quasiparticles are simply the allowed states on the annu- 
lus. We need only consider configurations with at most 
k curves extending from the quasiparticle to the outer 
boundary. The value of any physical wavefunctional on a 
configuration with more than k such curves can be related 
using Pfc+i to its value on a configuration with fewer than 
k + 1 such curves together with some curves which begin 
and end at the inner boundary of the annulus (and some 
which do the same at the outer boundary). The latter 
are edge excitations, which we ignore for now and discuss 
later. Thus, we can classify quasiparticles according to 
the number n of curves which terminate at them. In the 
level-fc theory n = 0, 1, . . . , k. 

For concreteness, we focus on the level k = 2 theory, 
but the extension to other k is straightforward. Either 
n = 0, 1, or 2 curves can terminate at the inner boundary 
of the annulus. Let us choose two preferred points on 
the inner boundary and two more on the outer boundary 
where curves can terminate. We furthermore specify that 
it is preferable for curves to terminate at one of these 
points; a curve will terminate at the other preferred point 
only when the first one is already taken. This is not the 
most natural boundary condition, but it is convenient 
for counting states and drawing pictures. We will not 
allow bigons, curves which have both of their endpoints 
at the same boundary. As promised earlier, a discussion 
of boundary conditions and edge excitations will follow 
in a later section. 

The allowed states on the annulus are depicted in fig- 
ure ED Note that there are three species of excitations 
(including the vacuum) in which there are no curves ter- 
minating at the quasiparticle, four species of excitations 
in which there is one curve terminating at the quasipar- 
ticle, and two species of excitation in which there are two 
curves terminating at the quasiparticle. 

Let us denote the three states in figure 118b as 
|0), |0 r ), |0 r 2). We denote the four states in figure IT%b as 
|1), I It), \1 t ~i), |1t 2 ), an d the two states in figure fHfo as 
|2), \2jj)- The spin and two-particle braiding eigenstates 
are linear combinations of these: 

spinO : V2\Q r ) ± |0 r a),2|Q) - |0 r =) 
spin ±i : \2)t*\2h) 

spin e(^) : ±e^ 8 \l)-e^\l T ) 

Te -«W8|l T _ 1 ) + |l Ta ) 
e = ±l (115) 

(N.B. By 'spin' s, we mean here that the state has eigen- 
value e 27 ™ s under a spatial rotation by 2ir. This is entirely 
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FIG. 18: The 9 quasiparticle species (including the vacuum, 
or 'trivial particle') of the k — 2 theory. See text for an 
explanation. 



distinct from the SU(2) x SU(2) representation which it 
might carry.) 

These quasiparticles can be associated with corre- 
sponding SU(2) representations. Let us assume that ev- 
ery particle carries two SU(2) quantum numbers (Ji, J2), 
with ji = 0, h, 1. Let us decompose the product of 
these quantum numbers: ji ® j% = \ji — J2I, \ ji — J2I + 
1, . . . , max ( ji + j 2 , 1). This is almost what we would ex- 
pect for the Lie group SU(2), except for the upper cutoff 
k/2 which, in this case, is 1. In the next section, we 
will see the origin of this upper cutoff in the language of 
Kac-Moody algebras. Now consider the 9 allowed com- 
binations (ii,j2) an d decompose their ji ® j 2 products 
into irreducible representation of total j. Four of these 
products are j = h. These correspond to the four states 
with a single curve terminating at the inner boundary, in 
figure H8b. There are also two j — quasiparticles and 
three j = 1 quasiparticles. These correspond, respec- 
tively, to the first two states in figure 118b and the two 
states in figure [T5b together with the third one in 118b. 

For general k, the SU(2)k theory is expected to only 
have the representations j = 0, |, . . . , § , and the (k + l) 2 
quasiparticle species can be organized in the manner just 
described for k = 2. in the next section, we will show how 
this can be facilitated with some results from conformal 
field theory. Incidentally, this is why there are only four 
states of two j — 1/2 quasiparticles in SU(2)k theory: 
pairwise, they must form j = or j = 1. 

Let us now consider the braiding statistics of these 
quasiparticles. If we have two quasiparticles, their braid- 
ing statistics is determined by the quasiparticle spins. 
Note that on a compact manifold, if there are only two 
quasiparticles, then they must fuse to something topolog- 
ically trivial. On the annulus, they could fuse to form a 
non-trivial quasiparticle because the outer boundary can 
compensate. For instance, suppose we have two n = 1 
quasiparticles. On the sphere, they must fuse to form 



26 



n = (in fact, they must fuse to form the trivial parti- 
cle); on the annulus, they can fuse to form n — 2 since the 
two curves can terminate at the outer boundary, which 
must also have n — 2. In either case, performing a braid 
cannot change n of the composite. In fact, we could braid 
one quasiparticle around another and then fuse them or 
simply fuse them and rotate the fused particle. This must 
lead to the same result, which is clearly just a phase in 
the second approach. 

Suppose we have four n — 1 quasiparticles in the k = 2 
theory. There are 5 possible states of 4 quasiparticles on 
the sphere. To sec this, note that fusing two of the quasi- 
particles can lead to one of the three n — quasiparticles 
or one of the two n = 2 quasiparticles. The other two 
quasiparticles must fuse to form the same excitation since 
the aggregate of all four quasiparticles must be topolog- 
ically trivial. These five states are depicted in figure IT§1 
(This is just one basis of the five states. One could eas- 
ily draw a different set of five pictures but they would 
be related to these using P 3 .) It is clear from the figure 
that taking particle 2 around particles 3 and 4 trans- 
forms the first state into the second. It is also clear that 
this does not commute with, say, exchanging particles 2 
and 3. Thus, these particles exhibit non-Abelian statis- 
tics. The k = 1 case is special because P2 allows us to 
collapse all of these states into one state. In this spe- 
cial case, the statistics is Abelian, since there is only one 
p-quasiparticle state, so it can at most acquire a phase. 
For higher k, Pk+i is simply not as restrictive. The braid- 
ing matrices for these theories can be computed simply 
by drawing pictures and applying the projector relation 
Pk+i- 




FIG. 19: A complete, linearly independent (but not orthogo- 
nal) set of 5 states of two n = 1 quasiparticles in the k — 2 
theory. 



VI. RELATION TO 2D CONFORMAL FIELD 
THEORY 

At first glance, 2D (or 1 + 1 — D) conformal field the- 
ories would seem to be ill-suited to describe topologi- 
cal phases of 2 + 1-dimensional physical systems. Obvi- 
ously, the dimensionality is wrong. Furthermore, confor- 



mal field theories describe critical points, at which the 
spectrum is gapless, not stable phases with a gap. How- 
ever, upon closer inspection, conformal field theories do, 
in fact, have many of the necessary ingredients for ex- 
otic braiding properties. Consider the holomorphic (or 
right-moving) part of a conformal field theory. Since the 
spectrum is gapless, the two-point correlation function of 
a field $ varies as a power law: ($(z)$(0)) ~ 1/ z 2h . If 
h is not a half-integer, then this correlation function is 
multi- valued. A phase is acquired as one field encircles 
another, analogous to what occurs when one particle is 
taken around another in 2 + 1 dimensions. While a chiral 
two-point function will suffer from at most a phase am- 
biguity on the plane or sphere, it can have more severe 
non-uniqueness on higher-genus manifolds as can the chi- 
ral part of a correlation function, of four or more fields 
(i.e. a conformal block) on the plane itself. 54 In general, 
there will be a vector space (which is finite-dimensional, 
in the case of rational conformal field theories) of confor- 
mal blocks. Taking one field around another will lead to 
monodromy matrices rotating the conformal blocks into 
each other. 

A second propitious feature of conformal field theories 
is the existence of an operator product expansion. When 
two fields are viewed from a distance much larger than 
their separation, they appear much as a single 'compos- 
ite' field would or, rather, as a linear superposition of the 
different possible 'composite' fields: 

(z) fa (w) = J2 c m (* - w) h *- hi - h ' cp k (w) (116) 

k 

so a number of different 4>kS can appear on the right- 
hand-side of this equation. This is analogous to the 'fu- 
sion' of two particles in 2 + 1 dimensions: when two par- 
ticles are brought close together, they can fuse to form a 
single composite particle. There are different possibilities 
for this particle, so the result of fusion is, in fact, a linear 
superposition of composite particles. 

This is a particularly convenient decomposition for dis- 
cussing braiding. When a counter-clockwise exchange of 
z with w is performed, the coefficient of 4>k changes by 
a phase factor e m ( hh ~ hi ~ h i' . In other words, the sum 
on the right-hand-side is over eigenvectors of a counter- 
clockwise exchange. If we consider the four-point correla- 
tion function, there will, in general, be several conformal 
blocks, which transform according to monodromy matri- 
ces as a result of exchange operations. These matrices 
have eigenvalues e K% ( hh ~ hi ~ h} \ 

The resemblance between conformal field theory and 
Chern-Simons theory is not coincidental, but reflects an 
underlying relationship between the theories. It is easi- 
est to see the connection between the two by considering 
Chern-Simons theory for some semi-simple Lie group, G, 
on the 2 + 1 — D manifold D 2 x R, where R is identi- 
fied with the time direction. In Coulomb gauge, a§ = 
(where a is a Lie algebra index), the Hamiltonian van- 
ishes and there is only the constraint /" = 0. In the 
functional integral formulation, this is expressed as fol- 
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lows. The field ag only appears in the action linearly, so 
the functional integral over may be performed, yield- 
ing a 5-function: 



Da, 5{m) lo'xR^ *^ 



where i, j — 1, 2 run over spatial indices. The constraint 
imposed by the (5-function can be solved by taking 



= diUU~ 



(118) 



where U is a single-valued function taking values in the 
Lie group. Substituting this into the right-hand-side of 
(|117f) . we find that the action which appears in the ex- 
ponent in the functional integral takes the form 



k 
477 



k 

4-7T 

k 

47T 



D 2 xR. 



Hi {d l UU- 1 d (d J UU- 1 )) 



D 2 xR 



D 2 xR 



127T 



D 2 xR 



^(d.UU^dodjUU- 1 ) + 
tr (d l UU- 1 d 3 Ud Q U- 1 ) 
d j ti(d i U- 1 d U) + 

tr (d l UU- 1 d J Ud Q U- 1 ) 

SixR 

e^hi (d^U V~ x dJJ U- 1 d x U U'^W) 



The Jacobian which comes from the ^-function <5(/° ) is 
cancelled by that associated with the change of integra- 
tion variable from Da to DU. In the final line, we have 
integrated by parts the first term so that the integral is 
over the boundary, which is parameterized by the coor- 
dinate x\. The second term has an extra factor of 1/3 
resulting from the more symmetrical form in which we 
have written it. Though it appears to be an integral over 
the 3D manifold, it only depends on the boundary values 
of U. This action is the chiral WZW action, as we will 
see momentarily. 

The equation of motion is simply 



d (U^diU) = 



(120) 



The equation of motion only deals with the restriction of 
U to the boundary, S 1 x R, which is one way of seeing 
that the action is independent of the continuation from 
S 1 x R to D 2 x R. This equation is a slightly obscured ver- 
sion of the chiral WZW equation of motion. The Chern- 
Simons action has a vanishing Hamiltonian which is why 
the current U~ 1 d\U is independent of time. With a more 
complicated boundary condition, we could impose non- 
trivial dynamics at the boundary, which would lead to 



a similar "off-diagonal" derivative structure in the equa- 
tion of motion, but with 8q — > dt ~ d\, d$ — > 8q + d±: 



d- (U^d+U) = 



(121) 



which is the usual equation of motion of the right-handed 
part of the WZW model. 

Either of these equations, (|12(J|) or i|121|) . states that 
the currents 



ja = tr ( T a U- 1 d 1 U) 



(122) 



are free, where the T a s are the generators of the Lie al- 
gebra in the adjoint representation. Consequently, they 
obey a Kac-Moody algebra, which allows us to alge- 
braically compute the braiding properties of primary 
fields. The OPE of the currents yields the Kac-Moody 
algebra: 



J a (xi) J b (0) 



k 5 



ah 



f abc J c {0) 



Xx 



(123) 



where f abc are the structure constants of the Lie algebra. 
The energy-momentum tensor of this theory is of the 
Sugawara form: 



T : 



1/2 
k + C A 



J a J a 



(124) 



where Ca is the quadratic Casimir in the adjoint repre- 
sentation if the highest root is normalized to length 1. 
A field which transforms in representation r of the 
group G and is primary under the Kac-Moody algebra, 



J a ( Xl )cp (r) (0) 



XI 



(125) 



has, according to (|124fl . dimension 

rta> *>«<«> = si^T 

_L_ C r 

x\k + C A 



(126) 



For the case of SU (2)^, this means that a spin j primary 
field has dimension hj = j(j + l)/(fc + 2). 

A restriction on the allowed j's in the SU(2)k theory 
can be found by expanding J a {x{) in modes 



J>i) 



s ja g— i(m+l)xi 



(127) 



Then the operators I a = Jq form an su(2) Lie alge- 
bra. Hence, 2Jg has integer eigenvalues in any finite- 
dimensional unitary representation. Similarly, I 1 = 
(Jl + Ji x ) /V2, P = (J? + J 2 _ x ) /V2, I 3 ee I k- J$/y/2 
also form an su(2) Lie algebra. Consider a spin j highest 
weight state \j, m = j), with I 3 \j, m = j) = j\j, m = j). 
Then 

< (j, rn = j\I + r \j, m = j) 
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(j,m = j\ 



I+J- 



\j, m = j) 



= (j, m=j\k- 2I 3 \j, m = j) 
= k-2j (128) 

Hence, the SU(2)k theory has particles which transform 
under the j = 1/2, 1, . . . , k/2 representations. These par- 
ticles have 'spins' hj = j(j + l)/(k + 2), which deter- 
mine how they transform under a spatial rotation. They 
have fusion rule ji <g> j% — ®ja with \j\ — ja| < ja < 
min(ji + j2, k — j\ — ja) (the derivation may be found 
inSi). Together, these completely determine the braid- 
ing properties of the particles. 

A spin-j primary field in the SU(2)k conformal field 
theory corresponds to a Wilson loop in the spin-j rep- 
resentation in the SU(2)k Chern-Simons theory. Thus, 
we can calculate the braiding properties of Wilson lines 
in Chern-Simons theory using the OPE in the associated 
conformal field theory. Consider the spin- 1/2 represen- 
tation of SU(2) k . Then, 



#/a (*) </ 2 (0) = * 



-3/(2fe+4) ,a/3. 



Cz 5/(2k+i) e «6 T . ^ (()) (129) 



where a, /?, 5 are spinor indices, rs are Pauli matrices, 
and C is a known constant. Thus, there is a 2 x 2 braid 
matrix associated with an exchange and its eigenvalues 
are - e 3*i/ t e ^i/(2k+i) _ Notej however, that there is 
an extra half- twist or a deficit of a half- twist, respectively, 
associated with each of these operations (i.e. the framing 
associated with them is different from the standard one 
which we chose earlier), so that the desired eigenvalues 
are actually - e ™/( 2fc + 4 ) , e 37ri/(2fe+4) _ 

This result can be used to calculate the value of a 
contractible, unknotted Wilson loop, according to the 
derivation of the previous section: 

d = 2 cos (--?—) (130) 



VII. EDGE EXCITATIONS 

In our previous discussion of the Hilbert space of dou- 
bled Chern-Simons theory in the presence of a mani- 
fold with boundaries or in the presence of quasiparti- 
cles, we explicitly forbid bigons with endpoints at the 
same boundary. Now, however, it is time to let bigons be 
bigons. We also remove the constraint which fixed the 
endpoints of curves to lie at marked points. As we will 
see, the boundary excitations can be understood in terms 
of the conformal field theories of the previous section. 

In general, the endpoint of a curve which terminates 
at a boundary will be described by a wavefunction ip{9). 
Previously, we required ip(0) = 5(8 — 9 ). This can be 
viewed as the correct ground state in the presence of a 
particular boundary Hamiltonian. The most important 



terms in the boundary Hamiltonian will be of the form 
H = Y J {PLI + V) (131) 

i 

In other words, there will be an energy penalty r\ for 
each endpoint and a kinetic energy proportional to the 
square of the angular momentum, Li = d/d6i. Thus, 
the eigenfunctions in the sector with n endpoints will be 
angular momentum eigenstates 



</>(0i,0 2 ,...,0„) 



(132) 



where 9±, 6%, . . . , 6 n are the angular positions of the n 
endpoints. These may be endpoints of bigons or of curves 
which ultimately terminate at other boundaries. 

The allowed values mi are determined by consistency 
with the surgery relation of the theory. In general, they 
are not integers. Let us consider, as examples, k = 1,2, 
corresponding to d = —1, \[2. First, consider states with 
a single endpoint in the d = —1 theory. Using P%, we 
see that the wavefunction = e lm0 must satisfy the 
relation ip(6) = —ip(0 + 4tv). Hence, m € Z ± 4. States 
with more than one endpoint are equivalent to states 
with either zero or one endpoint and some number of 
bigons. Under a 2tt rotation, the bigons are unchanged, 
so the total angular momentum M = X)i m i must satisfy 
tog Z or m 6 Z ± 4, respectively, for an even or odd 
number of endpoints. 

Now, consider the k — 2 theory. For a single endpoint, 
P 3 dictates that ip(9 + 2tt) = e 2ms ip(9 + 2tt), where s 
takes one of the four values s = ± (^p) , as in eq. (|115|l . 
In other words, m G Z + s. For two endpoints, it dictates 
that s = ±i. States with more than two endpoints can 
always be related using P3 to a state with 0,1, or 2 end- 
points and some number of bigons. The bigons are, of 
course, invariant under 9 — * 9 + 2-7T, so the constraint on 
the angular momentum is the same as above, depending 
on the number of endpoints modulo 3. 

Thus, we see that the boundary states of the doubled 
Chern-Simons theories can be understood as massive de- 
formations of the corresponding achiral conformal field 
theories which we discussed in the previous section. If 
we were dealing with undoubled chiral topological theo- 
ries, the corresponding edge excitations would be chiral 
and, therefore, necessarily gapless. Achiral theories must 
be tuned - so that 77 = in (|131f) - in order to be gap- 
less. The basic state counting is the same, however. In 
the SU(2)x x SU(2) 1 case, we have the spin sector - or 
states with an even number of endpoints - corresponding 
to the towers of states: 



Ji 



rL 
— 711 

L 

— Urn 



■ Ji 



J\...J\\Q) 



(133) 



J\...j\^ /2 ^n (134) 



The rii , ^ correspond to the angular momenta differences 
between the two endpoints of the associated bigons. <j>y 2 , 
4>i/2 are the SU (2) doublet primary fields. We also have 



the spin ± j sectors - states with an odd number of end- 
points - corresponding to the towers of states: 

jf ni ...jf nm jv..jVf /a |o) ( 135 ) 

jf ni ...jf„ m A...Jf,,^ /2 |0) (136) 

The fields </>fy 2 ' ^1/2 have spins ±^, confirming the cor- 
respondence. 

SU(2) 2 x SU {2)2 has the states with zero endpoints 
modulo 3: 

jf m ...jf„ m J\...J\\Q) (137) 
jf ni . . . Jf „ m J\ . . . J\ 0? /a # /a |O) (138) 
jf ni ...jf„ m jf^-J^f^O) (139) 

The states with one endpoint modulo 3: 

J-n,---J- nm A--- A^/al°) ( 141 ) 
J-m ■ ■ ■ A m A • ■ ■ A^^/al°) ( 142 ) 
Jf Bl . . . Jf „ m J\ . . . J R lq 0f 0[ /2 10) (143) 

(144) 

have spins with are an integer plus ±3/16, ±1/16 since 
<f>y2 nas s P m ±3/16 and (/) R ' L has spin ±1/2. The states 
with two endpoints modulo 3: 

Jf ni ...jf„ m J R ll ...J R lq ^\0) (145) 
Jf ni ...jf nm J\...J R lq 4>m (146) 

have spins ±1/2. 

Continuing in this way, we could construct the edge ex- 
citations for any of the SU(2)k x SU{2)k theories: they 
are simply massive theories corresponding to the associ- 
ated achiral conformal field theories. The allowed weights 
of the primary fields follow from an application of the 
projector relation Pk+i, as we have shown explicitly for 
the cases of k = 1, 2. 



VIII. TOWARDS MICROSCOPIC MODEL 
HAMILTONIANS 

Our formulation of doubled Chern-Simons theories in 
terms of Hilbert spaces spanned by configurations of 
multi-curves on surfaces is compact and elegant. How- 
ever, its virtues are not purely esthetic, but also include 
its natural connection to microscopic models which give 
rise to these phases. In section ITTT1 we showed how they 
could arise from mean-field solutions of a variety of micro- 
scopic models of interacting electrons, but did not show 
that these were the true ground states of any particular 
Hamiltonians. Here, we take a different tack and take 
some steps towards a more direct connection between 
microscopic models and the loop space formulation of 
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FIG. 20: S z basis states of s = 1/2 spins on the triangular 
lattice can be represented in terms of loops on the honeycomb 
lattice. 

doubled Chern-Simons theories. As we observe in this 
section, many systems admit a loop space description. 

Consider a system of s = 1/2 spins on a triangular lat- 
tice. Let us work in the S z basis in which every spin takes 
a definite value |, |. Let us represent these basis states in 
terms of the domain walls which separate clusters of up- 
and down-spins. For every configuration of domain walls, 
there are two spin configurations which are related by a 
reversal of all spins. The domain walls lie on the links 
of the dual honeycomb lattice and the spins sit at the 
centers of the faces of the honeycomb lattice. It is clear 
that these domain walls can neither terminate nor cross. 
Thus the Hilbert space of a triangular lattice spin system 
is of precisely the desired form, as depicted in figure |2*U1 
However, only very special Hamiltonians will lead to a 
ground state which obeys relations such as (|114fl . 

For the initial members of our sequence of theories, 
d = ±1, such Hamiltonians can be written in a simple 
form. Consider theory B, the d = 1 theory. The Hamil- 
tonian H = h^^f requires every spin to point in the 
x-direction in spin space or, in other words, in an equal 
linear superposition of the states | |) and | J,). Flipping a 
spin causes the dark and light bonds of the corresponding 
hexagonal plaquette to be exchanged in figure |2*U1 This 
can either create a new contractible loop, erase a loop 
of minimal size, deform a loop, or perform the surgery 
operation )( — > ~. 

Thus, we have found that a system of spins in a mag- 
netic transverse field is trivially equivalent to the sim- 
plest of our topological theories. This is a special case 
( J = = K = 00) of the duality between the transverse 
field Ising model 

h = -jJ2 s * s j + h Y. s ? ( 147 ) 

<i,j> i 

and the Z2 gauge theorj^&. 
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However, this is only the 'even' part of the theory. Do- 
main walls cannot terminate in the bulk of the system 
and they can only terminate in pairs at a circular bound- 
ary. Thus, the above Hamiltonian does not admit ex- 
citations such as the one depicted in figure [5] However, 
the following simple model due to KitaeA*2I contains both 
the even and odd parts of the theory. There is a spin- 
1/2 degree of freedom a z — ±1 on each link of a lattice. 
The lattice is arbitrary, but let's consider a honeycomb 
lattice, for the sake of concreteness. 



H = J^A, - J 2 J2 f p 



(148) 



where 



A, 



These operators all commute, 



[F p ,F p ,) = [A i ,A j ] = [F p ,A j ] = 



(149) 



(150) 



so the model can be solved exactly by diagonalizing each 
term in the Hamiltonian: the ground state |0) satisfies 
Ai\0) = -|0),F P |0) = |0). If we represent a z = 1 by 
colored bonds and a z = — 1 by uncolored bonds, then 
Ai\0) = — 10) requires chains of bonds to never end, 
while F P \Q) = |0) requires the ground state to contain 
an equal superposition of any configuration with one ob- 
tained from it by creating a new contractible loop, erasing 
a loop of minimal size, isotopically deforming a loop, or 
performing the surgery operation )(— > This is clearly 
the same as the transverse field spin model above, except 
that curves can now terminate, albeit with an energy cost 
2Ji. 

In a similar way, we can formulate a model which gives 
rise to the d = — 1 theory. 



H 



Jd, s y^ J F°' d 

p 



(151) 



where 



6 



F: 



A v = ^-aeU{v)^a 

F; ee at at at at at + h.c. 

+ vt< 7 2< J t (T i< 7 5_ rT >_ 
+ °l°2 CT 3 G A °t CT 6 

+ cyclic permutations 

.d x x x x „x „x T^i 



h.c. 
h.c. 
h.c. 



(152) 



The first term in the Hamiltonian is the same as in the 
d = 1 theory: it selects a low-energy subspace in which 
chains of up-spins never terminate. Hence, the config- 
urations of this low-energy subspace can be represented 
as closed multi-curves. In the d = 1 theory, isotopy, the 
condition d = 1, and the surgery relation can all be imple- 
mented with a single plaquette term. For d = — 1, how- 
ever, the ground state must contain an equal superpo- 
sition of isotopic configurations but superpositions with 



opposite signs of configurations which are related through 
surgery or the addition of a contractible loop. These con- 
ditions are enforced by the second and third terms in the 
Hamiltonian. (The fourth term in F* is actually the re- 
sult of two surgeries and, therefore, carries a factor of 
d 2 = 1, the same as the isotopy moves with which we 
have grouped it.) The vertex and plaquette terms in the 
Hamiltonian commute with themselves and each other. 
Hence, this Hamiltonian is exactly soluble. Following 
steps similar to those of the previous paragraph, we see 
that this Hamiltonian implements the d = — 1 theory and 
supports semionic excitations. 

We expect that other physical systems can give rise to 
this type of structure and, with some luck, topological 
phases corresponding to doubled Chern-Simons theories. 
Dimcr models have a natural representation in terms of 
multi-curves. In these models, it is assumed that that 
there are spins located at the sites of a lattice and that 
each spin forms a singlet dimer with one of its nearest 
neighbors. Thus, there are dimers on the links of the 
lattice which satisfy the following condition: there is one 
and only one dimer touching each site. A dimer covering 
is not, of course, composed of closed curves, but the tran- 
sition graph between two dimer covering is: one considers 
some fixed reference dimer covering R and superposes it 
on the dimer covering of interest C. Where R and C 
coincide, we erase the dimers (or, perhaps, think of them 
as a minimum size closed loop). The remaining dimers 
form closed loops, with dimers from R and C alternating 
as one travels along the loop, as shown in figure |^ A 
dynamics is now needed which assigns a weight d to con- 
tractible loops, allows the loops to deform isotopically, 
and enacts Pk+i- The Kivelson-Rokhsar 21 Hamiltonian 
on the triangular lattice^ does this for the d = 1, k = 1 
theory. Some ideas about implementing the higher-A: the- 
ories in dimer models are discussed in£&. 

Josephson junction models^Si admit a similar de- 
scription. In these models, superconducting islands 
(which, in the simplest incarnation, live on the links of a 
lattice) are connected by tunneling junctions. The charge 
is effectively restricted to take the values 0,1, which cor- 
respond to the presence or absence of a curve on that 
link. The Hamiltonian further requires these curves to 
be connected and non-terminating. More elaborate mod- 
els allow for charges 1, 2, . . . , N, which can be mapped to 
labeled curves. These may be useful for the implementa- 
tion of k > 1 theories since a curve carrying the label n 
might be a way of representing n curves. 

Of course, finding a model with a loop gas representa- 
tion is only the first step. The model must also (a) assign 
a value of d to each closed contractible loop and (b) en- 
force the associated Jones- Wenzl projector relation. It 
might seem improbable that any realistic Hamiltonian 
would impose precisely the right value of d, much less 
impose the corresponding Jones- Wenzl projector, which 
can be quite complicated, as may be seen from figure 
ITfil However, if a model incorporates the correct value of 
d = 2cos(7r/fc + 2), the uniqueness of the corresponding 
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FIG. 21: By superimposing a dimer configuration (dark lines) 
on a reference dimer configurations (dotted lines), it can be 
represented by the depicted multi-loop which consists of al- 
ternating solid and dotted lines. 



Jones- Wenzl projector implies that a generic perturba- 
tion stands a good chanc e of driv ing the system into the 
corresponding SU(2)k x SU(2)k topological phase. 



IX. A 'PLASMA ANALOGY' USING LOOP 
GASES 

Some insight into the physics of the wavefunctions 
discussed above, their representation of the topological 
structures discussed in this paper, and the difficulties 
associated with imposing the Jones- Wenzl projector on 
them with a local perturbation can be gained by con- 
structing a "plasma analogy" for the ground state wave- 
functions of our topological field theories. As in the case 
of Laughlin's plasma analogy for Abelian quantum Hall 
states, the idea is to map the squared norm of the quan- 
tum ground state in a 2 + 1-dimensional system to the 
partition function of a classical 2-dimensional system. 
Equal-time correlation functions in the former will be 
correlation functions in the latter. However, unlike in 
Laughlin's plasma analogy, in which the squared norms 
of ground states were mapped to the partition functions 
of plasmas, we map them to loop gases. In order to effect 
this mapping, we will need to consider classical models on 
lattices. Such a short-distance regularization may seem 
unnatural and unnecessary from the point of view of the 
preceding discussion of doubled Chern-Simons theories. 
However, it will prove to be very natural when we turn 
to 'microscopic' lattice models which support topological 
phases. 

Let us consider the unnormalized ground state on the 



sphere for a given value of d. It is given by 

|0>=X>»«|{7}) (153) 
{7} 

where n c (pf) is the number of loops in the configuration 
7. Its norm is 

(0|0)=]Td 2 "^>({ 7 }|{ 7 }) (154) 
{7} 

In other words, it is a sum over all possible loop config- 
urations, weighted by a loop fugacity d 2 . 

Two types of classical statistical mechanical models 
have very similar sums appearing in their partition func- 
tions, the ferromagnetic g-state Potts model and the 
0(n) model. The details of these mappings can be found 
in Ref. [57] and (the former) in Ref. l58l nevertheless we 
shall briefly review them here. 

We consider first a set of 0(n) models which give rise 
to loop gases 5 «. These models have partition function 




(155) 



where the S{S are ./V-dimensional unit vectors, Kn is the 
surface area of the iV-sphere, and 

~[3H = J2H l + xSi -Sj) (156) 

or, simply, 

e -W = YKl + xSi-Sj) (157) 

{i,3) 

While the Hamiltonian can take unphysical imaginary 
values for x > 1, the model is still perfectly well-defined. 
Let us work with this model defined on the honeycomb 
lattice. 

The high-temperature series expansion of this model 
is obtained by expanding the product i|157|) . Consider 
any given term in this expansion. For each link on the 
lattice, there will either be the factor 1 or xSi ■ Sj. If it's 
the latter, then this term will vanish upon integration 
over Si and Sj unless there is another factor in the term 
which contains Si and a factor containing Sj. Suppose 
the factor fulfilling the former requirement is xSi ■ Sk- 
Then a factor containing §k must also be present. Con- 
tinuing in this way, we see that if there are any bonds 
of the lattice for which factors of the form xSi ■ Sj ap- 
pear rather than 1, then these bonds must form closed 
loops or else the term will vanish. Clearly, there is no re- 
quirement of close-packing. There is one non-vanishing 
term in the expansion in which there are no bonds. By 
choosing the honeycomb lattice, we have ensured that 
the loops cannot cross. 

For each vertex on such a loop, we have a factor of 

yg 5?5 f=i^ (158) 
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Thus, for each closed loop, we obtain a factor of 

^ N N "' N 

-(£)'" < 159 > 

Hence, the partition function is 

^ = E(|)V (160) 

G 

where b is the number of bonds and I is the number of 
loops. In this way, for x — N and N — d 2 , we have a 
statistical mechanical model whose partition function is 
the squared norm of the ground state of SU (2)k x SU(2)k 
Chern-Simons Gauge Theory. 

The g-state Potts models also have loop gas represen- 
tations, but, as we will see, the loops are fully-packed in 
this case. With no wiggle room, isotopy is impossible. 
It is conceivable that this makes life easier since there is 
no need to impose isotopy invariance as a condition on 
low-energy states. Thus, it may prove useful to consider 
microscopic models of this form. 

The Hamiltonian of the ferromagnetic (/-state Potts 
model is given by 

where (Tj = 0, 1, . . . q — 1, J > 0. With the help of the 
identity exp (Jc^Oj) = 1 + vb aiaj where v = exp(J) — 1, 
the partition function for for this model can be written 
as follows: 

z Potts = j2 e ~ 0H = e n i 1 + vs °^ ) ■ ( i62 ) 

M M (hi) 

Expanding the product in (|162|) can be interpreted graph- 
ically: every time v8 airj - is chosen for a pair of neighbor- 
ing sites i and j, the corresponding bond is colored; the 
choice of 1 results in an empty bond. Due to Kronecker 
(5-symbols, all sites belonging to the same cluster must 
have identical values of spins a. Summing over all possi- 
ble spin configurations {a} we then obtain 

Zpotts = 5>Y (163) 
Q 

where b is the total number of occupied bonds and c is the 
number of clusters (including isolated sites) . The sum is 
now performed over all configurations Q of such clusters. 
This is a so-called Fortuin-Kasteleyn or random cluster 
representation^. 

A so-called polygon decomposition^ lets us relate this 
to a loop gas on the surrounding lattice (vertices of 
the surrounding lattice are the midpoints of the origi- 
nal bonds). We then think of an occupied bond as a 
double-sided mirror placed at the site of the surrounding 



lattice. If a bond is not occupied, then its dual bond is 
considered a mirror. Thus every site of the surrounding 
lattice gets one of the two possible mirrors. We then use 
these mirrors to construct paths as shown in figure 1221 
Since these paths have no sources or sinks, they always 
form loops that either surround the clusters or are con- 
tained inside clusters (in the latter case, the loops can be 
thought of as surrounding dual clusters) . The number of 
loops I is given by £ — c + / where / is the number of 
faces, i.e. the minimum number of occupied bonds which 
have to be cut in order to make each cluster tree-like 
(essentially the number of "voids" which are completely 
contained within clusters). If we use the Euler relation, 

b + c- f = const. (164) 

then I = 2c + b, or c = (£ — b)/2. Hence, we can rewrite 

^Potts as 

z Potts = 5>y £ - b)/2 

- < i65) 

If v = ^/q - i.e. if the Potts model is at its self-dual point, 
then 

«; Dual = E(v^ ( 166 ) 

G 

This appears to be a lattice regularization of the norm of 
the ground state wavefunction above if d 2 = y/q, apart 
from the full-packing condition: that is to say that every 
bond of the surrounding lattice belongs to a loop, no 
bonds are left empty, as seen in figure [221 Since there 
is no corresponding notion in the continuum, it is not 
completely clear what the connection is to our earlier 
discussion of the ground states of doubled Chern-Simons 
theories. 

The O(N) and self-dual Potts models are quite differ- 
ent, but they share one curious feature, which is of in- 
terest to us: for d > 2, which corresponds to N > 2 and 
q > 4, respectively, all correlation functions are short- 
ranged. The underlying physical reasons are unrelated: 
the 0(N > 2) models do not have an ordered phase by 
the Mermin-Wagner theorem, so their correlations are 
short-ranged. In the ordered phase, loops become long 
and wander around the system, as happens for N = 1, 
i = l, which is deep within the ordered state of the 
Ising model. For N = 2, there is an algebraically or- 
dered phase, which includes x — 2, so loops are still able 
to wander about the system. The q > 4 self-dual Potts 
models, on the other hand, are short-ranged because they 
have a first-order phase transition. The ordered and dis- 
ordered phases have only small loops, and loops can grow 
large and wander only at a second-order phase transition 
point, as occurs for q < 4. For q > 4, the phase transition 
is first order, and loops never grow large§£. 

The basic phenomenon is that for d large, it is highly 
favorable to have as many contractible loops as possible. 
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FIG. 22: A typical cluster configuration for the Potts model 
is shown by dashed lines. Spins belonging to the same clus- 
ter take the same value, which must be summed over the q 
possible values, as described in the text. Clusters can be rep- 
resented by loops on the surrounding lattice, shown by solid 
lines. 

On the lattice, the proliferation of small loops takes up 
all of the space, leaving no room for loops to stretch 
and fluctuate. There is no corresponding effect in the 
continuum, where there is no shortest length. 

This makes it difficult for any local term on the lattice 
to enact the Jones- Wenzl relation on essential loops, i.e. 
long loops which encircle non-trivial topological features 
such as handles or quasiparticles or terminate at quasi- 
particles. Such loops have very low probability of ever 
approaching each other. In order to circumvent this dif- 
ficulty, we need to specify the short-distance dynamics of 
a candidate lattice model so that it more faithfully rep- 
resents the topological relations of the continuum theory. 

X. LATTICE MODELS, CHIRAL ANOMALIES, 
AND DOUBLED THEORIES 

From the preceding sections, we see that there is poten- 
tially a natural relation between doubled Chern-Simons 
theories and lattice models whose configurations can be 
represented in terms of loops. In this section, we note, 
on general grounds, that lattice models should yield dou- 
bled, rather than undoubled theories. 

Consider a one-dimensional electron system at finite 
density. Such a system is a metal. In the continuum, the 
numbers of left- and right- moving electrons are indepen- 
dently conserved at low energies. Their sum is the total 
charge, p, which must be conserved under all circum- 
stance. Their difference is the current, j, which suffers 
from the chiral anomaly, 

d»fn = E (167) 

where jp — [j, p) and E is the electric field. Since an 
electric field generates a current, the numbers of left- and 
right- moving electrons are not independently conserved 



in the presence of an electric field. In particular, in a 
homogeneous state, 

% = E (168) 

In other words, the chiral anomaly is simply the state- 
ment that a metal conducts electricity. To be completely 
precise, it is more than that, since it implies perfect con- 
ductivity, but only a little more since in a translationally- 
invariant system in which all particles have the same 
charge, a metal must be a perfect metal. 

In a crystal, which has only discrete translational sym- 
metry, this equation must be modified by the formation 
of bands. Instead, we have, 



v(k) = V fe e(fc) (169) 

with j given by a summation of v{k) over occupied states. 
Since e(fc) is a periodic function of fc, the current can- 
not increase linearly with E. Averaged over (sufficiently 
long) time, it is, in fact, constant. Said differently, if we 
follow the motion of the spectrum, then as right-moving 
states flow off to the right, just as many left-moving states 
flow in from the left because the number of these states is 
conserved and there is no net spectral flow. However, in 
the continuum, there is an infinite Dirac sea, and right- 
moving states can flow to the right (and be replenished 
by the infinite Dirac sea) without a compensating flow of 
left-moving states. 

This is not a great concern in the case of real met- 
als, which are not perfect metals because they are not 
translationally-invariant. They have finite conductivity 
as a result of scattering by impurities, phonons, etc., so 
they could never satisfy (|167[) anyway. However, there 
is a serious problem in chiral systems, since they must 
truly be perfect metals. Chirality implies that their cur- 
rents can not be degraded by backscattering. However, 
we have just seen that it is not possible to have a perfect 
metal, in the sense of <|167ll . on the lattice. This implies 
that a system on a lattice cannot be a chiral metal. 

Since there is no chirality in 2 + 1-dimensions 69 , these 
arguments do not apply directly. However, according to 
the arguments of section IVT1 the ground state wavefunc- 
tions of the TQFTs discussed here are related to cor- 
relation functions in associated conformal field theories. 
For chiral TQFTs, these conformal field theories are chi- 
ral, and have chiral anomalies. However, such a theory 
cannot arise from a lattice model. Thus, we do not ex- 
pect chiral TQFTs to arise from lattice models. Doubled 
theories are more natural. 

This does not mean that chiral theories are impossible 
on the lattice, just that some way of evading the above 
logic is needed. In order to have a finite chiral anomaly, 
an infinite Dirac sea is needed. This is not so artificial in a 
2 + 1 — D model since the excited states (which go beyond 
the mapping of the ground state to a 2D theory) can serve 
as the necessary reservoir. Such degrees of freedom would 
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seem to be missing from the pre-Hilbert spaces which 
we have been using in this paper, but a more general 
class of models might have the requisite structure. One 
simple possibility, motivated by the domain- wall fermion 
proposal for chiral lattice fermionsffi, is if the system has 
a finite extent in a third spatial direction (which one is 
free to envision as an internal degree of freedom) . 

XI. DISCUSSION 

Our purpose in this paper has been to explore the 
Hilbert spaces of a set of P, T-invariant topological field 
theories. These Hilbert spaces are most neatly under- 
stood in terms of a combinatorial construction which al- 
lows the distinct Hilbert spaces of the theories on differ- 
ent manifolds and with different quasiparticle numbers 
to be discussed in a unified formalism. The power of this 
construction derives from its reduction of the structure 
of these Hilbert spaces to a set of local rules for multi- 
curves on surfaces. Aside from its mathematical beauty, 
this is a Good Thing since it gives us some clues about 
what types of physical systems - which should, after all, 
be described by local rules - can manifest phases which 
are described by these topological field theories. The two 
most salient features are that the system should admit 
a description in terms of configurations of multi-curves - 
e.g. domain walls, chains of spins, etc. - and that higher- 
level theories should have longer (but still finite) ranged 
interactions or larger building blocks (e.g. higher spins). 

The essence of these phases is their ability to support 
excitations with non-trivial braiding statistics, which in 



almost all cases is non-Abelian. This property makes 
these phases relevant as a setting for topological quantum 
computation2L2&. It also makes them difficult to detect 
experimentally. Ideally, one would like to create excita- 
tions and manipulate them in order to perform braiding 
operations and Aharonov-Bohm interferometry. A more 
indirect way may be through the observation of a bro- 
ken symmetry phase which may be in close proximity 
to the topological phase of interest. If a finite density 
of Abelian anyonic excitations is created (by doping an 
insulator, say) then they will superconduct in zero mag- 
netic fielc&S. If the same is true with a system of non- 
Abelian anyonic excitations (see ref. for an argument 
that it is), then the observation of such an exotic super- 
conducting state may reveal the existence of a nearby 
topological phase. (Depending on one's perspective, one 
might argue that one or the other is the more interesting 
phase.) 
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